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é) “Paul, Why are you doing this special session ” & 


¢ Multi-purpose: ee 
— Educate those new to field of laser . 


based measurements (students) Hypersonic Nonequilibrium 
— Let collaborators from other fields le — _ 
know what is possible to measure 
— Provide some updates and ideas Edited by 
for those already working in 
measurement technology 


¢ Special session is coordinated 
with the launch of a new AIAA 


Progress Series Book 
— Pedagogical treatment : 
— Came out of a Fluids TC Working Ya eee 

group led by E. Josyula r | 
— Our chapter: Molecular-Based 

Optical Diagnostics for Hypersonic 

Nonequilibrium Flows ; Z 


* 


Timothy C. Lie 
. Volume 247 


uwen, Editor-in-Chief 
ume 2 
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About the Book 


The high-temperature environment poses an unusual challenge in understanding the basic physics of hypersonic flight. 
A lack of such understanding can lead to risks and uncertainties in the design of aerospace vehicles. Hypersonic 
— . Nonequilibrium Flows documents recent, unprecedented scientific advances in the field of nonequilibrium processes 


Most Recent Book 


? Available Books nat for aerospace applications. These advances have been driven primarily by interest in space access and exploration, or 
in developing military technologies involving hypersonic flight regimes. 
In the modeling of hypersonic flows, the last decade has witnessed a reexamination of fundamental principles of 
PRICES kinetic theory and quantum chemistry to describe the kinetic and thermal states, respectively, of the individual gas 
particles. Modern aerospace programs where nonequilibrium energy transfer processes play a major role may broadly 
PDF be categorized as exo- and endo-atmospheric. The exo-atmospheric flight programs consist of Earth and planetary 
Member: $79.95 reentry programs, as well as access to space programs for applications that include space exploration, flight 
List: $116.95 experiments and demonstrations, missile defense, and the nascent space tourism industry. The endo-atmospheric 
| Add to Cart | flight programs are primarily motivated by hypersonic military applications requiring high-precision engagement for 
SS tactical, theater and strategic defense, as well as applications involving intelligence, surveillance, and reconnaissance. 
Hardback The nonequilibrium processes considered in this volume are generally associated with flight Mach numbers between 7 
Member: $89.95 and 25, where the shock-layer temperatures range from 3000 to 25,000 K. 
List: $129.95 
| Add to Cart | Hypersonic Nonequilibrium Flows includes fundamental governing equations of nonequilibrium fluid transport and 


computational approach to calculation of rates and cross-sections in quantized energy states; DSMC approach; 
radiative heat transfer; a CFD perspective; surface chemistry, with additional chapters on high enthalpy facilities and 
the associated diagnostic techniques. 
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| Outline for 1° Talk S& 


Introduction and Application Considerations 
for Optical Diagnostics in Hypersonic 
Nonequilibrium Flows 


¢ Temperature, equilibrium and nonequilibrium 

¢ Characteristics of hypersonic nonequilibrium flow 
— What needs measuring? 

¢ Advantages of spectroscopic measurement 

¢ Aspects of a good measurement technique 

¢ Choosing which measurement technique to use. 


é Temperature and Equilibrium 


¢ Temperature 


— Classical thermodynamics 
¢ Heat flowing from object Ato B DT, >T, 
— Statistical thermodynamics 
¢ Describe by Boltzmann distribution 
¢ e.g. Maxwellian Velocity Distribution 
width ~ VT 


http://en.wikipedia.org/wiki/Temperature 


@ Equilibrium Greg L at the English language Wikipedia 
— Classical thermodynamics 
¢ Thermal equilibrium: T, = T, 
— Statistical thermodynamics 
¢ One temperature can describe all the 
different energy modes (via Boltzmann) 
— Chemical equilibrium: 
¢ Reactants and product concentrations 
don’t further change with time x 


population 


Ta> Tp 


Translation 


Rotation 


Vibration 


Electronic 


Energy Modes 


All these 
states are 
quantized! 
(ee ee 


\ \ 
Chm 
\ / j 


= — 


These can all be described by a temperature: 


—Trrang 


trans? ° rov’ 


Tip I, 


elec 


— Thermal equilibrium: these T’s are all equal 


S 


| Nonequilibrium S& 


¢ Thermal Nonequilibrium: 
— Different energy modes described by 


different T’s Rotation 

© Cee eee Population 
a nae , ; Distributi 

— Or, the distribution of population can’t eure 
be described by Maxwell-Boltzmann 


Statistics (“non-thermal” distribution) | 
| he 


¢ Chemical Nonequilibrium: 
— Products and reactants are in a state 
Rotational Quantum Number, J 


of change Computed using LIFBASE 
— If held at those conditions (P T)the = “wae Set imemational 


composition will change. 
¢ Equilibrium occurs via collisions 


é Hypersonic Flows and Nonegq. S& 


¢ Hypersonic vehicles fly at high altitude (low P) 
and high speed (M > 5) resulting in gas and 
vehicle heating (high 7) 

— Collisions rate scales as PVT 
— Low collision rate causes nonequilbrium 
— At STP for O,: Trans: 10; Rot: 10; Vib: 20,000; 
Dissoc: 200,000 collisions to equilibrate! 
° 10 coll.~1ns; 20,000 coll. ~2 us at STP 
¢ Hypersonic flows are also characterized by: 


— Steep gradients (shocks, boundary layers, 
expansions) 


— Combustion (in scramjets) 
— Dissociation / recombination 


é Some Hypersonic Nonequilibrium Flows S& 


Vine T, 


trans rot 
= am —— 


e Shock wave: 


thermal 


equilibrium ‘harmal 
non-equilibrium 


trans ~ ‘rot = 


X 
¢ Expansion: | thermal ee 
p chemical High PT. chemical Low al NO 
equilibrium non-equilibrium High V,M 


¢ We can measure, quantify these! 


| Parameters to measure S& 


¢ Typical parameters of interest in Hypersonic 
Nonequilibrium Flows: 


— Velocity, Pressure, Density, Temperatures, Species 
¢ Non-thermal population distributions 
— How is energy in flow partitioned in terms of: 
¢ thermal energy (T) 
¢ chemical energy (species) and 
¢ kinetic energy (V)? 


é Some advantages of Spectroscopic aa 
Measurement Techniques 
¢ Non-intrusive 
— No probes, particles perturb the flow 
¢ Remote detection 
— Can be used in hostile environments 
Can measure multiple parameters, states: | q 
— Sometimes simultaneously 
— Can be instantaneous (<10 ns) and high freq. (MHz) 
— Can be precise, accurate, high spatial res. (2D, 3D) 
¢ But... No one measurement technique can meet all requirements 
— Can have more complicated theory, expensive and 
complex equipment, more difficult data analysis 
— May require seeding of (sometimes toxic) gases, 
more expertise, slower data acquisition, analysis 
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é Scope of the Content Ss 


¢ Content focuses on molecular based, quantitative, 
mostly spectroscopic measurement techniques 
previously applied in or may be applicable to non- 
equilibrium flows. Many examples in equilibrium. 

— Content of lectures and manuscript are not exhaustive. Not 
all the ‘first’ or “most important’ references are given. Instead 
they there is an emphasis on describing how the techniques 
work and how they can be used to measure. 

— Also there is an admitted bias towards the authors’ own work. 


¢ Content excludes the many relevant particle-based (PIV, 
LDV), probe-based and surface-based measurement 
techniques that could be used to study non-equilibrium 
flows; mostly no ‘flow vis’. (No schlieren.) 


| What makes a good measurement technique? Ss 


¢ Areliable measurement system 
matched with well-understood physics 
to make a quantified measurement 
that meets a customer’s requirements. 


Explore each of these... 


o Typical Measurement System @y 


processing 


¢ Need to understand equipment and physics of 
process to make a measurement 
— Emission, Abs., LIF, Rayleigh, Raman, CARS, etc. 16 


| Measurement System: Lasers @Y 


e Lasers: pulsed vs. continuous wave (cw) 
— Continuous wave: lower power (mW to Watt) 


— Pulsed: 

10 ns, picosecond, femtosecond all commercially 
available 

‘Freeze’ the flow 

Much higher power vs cw (MWatt-GWatt, even TWatt) 
Frequency conversion (useful in linear techniques such 
as Raman, Rayleigh, LIF) 

Provide high powers for nonlinear techniques (CARS) 
Repetition rate: 10-100 Hz, 16 kHz, 1 MHz, 80 MHz... 


¢ What colors of laser light do we need? 


Colors of Light: Soectroscopy Ss 


Molecular response: vibrational 


electronic rotational 
sf i =< Increasing Frequency (Vv) 
. 10" v(Hz) 


1074 10°? 10°" 10'* 19'® 1o'4 10! 19!” 10° 10° 104 10 
| I | 1 | I I | | 1 | 
y rays X rays UV IR Microwave |FM AM Long radio waves 
Radio waves 
| 


1 | I | I th 18 | I 1 | I 1 

io io 10° 10° io* |} 0° io 10* 10° 107 10° 10° 10° 2 (m) 
' i 

Increasing Wavelength (1) = 


-_-" 


Visible spectrum Gace 
(from Wikipedia, “frequency”) 


. http://creativecommons.org 
Diode /licenses/by-sa/3.0/deed.en 


laser 
absorption 


LIF 


Raman/ 
Rayleigh 


500 600 700 


400 


Increasing Wavelength (A) in nm — 


¢ Many different colors can be useful for spectroscopy 


e As energy increases, more modes can be excited 
— LIF: can see electronic, vibrational and rotational structure 
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| Measurement System 2 


¢ Detectors: 
— Single point detectors: as 
Photomultipler tubes (PMT) and photodiodes 


— Camera technology is rapidly evolving: 


¢ CCDs provide ‘scientific grade’ (low noise, linear, etc.) 
¢ CMOS cameras are faster 


¢ Commercially available cameras at up to 10’s, 100’s of 
kHz steady streaming or MHz for short bursts or few pix 


— Fiber optics and fiber-optic bundles often used for 
either transmitting laser or collecting signal 
¢ Processing: 


— Some data processing codes are available (CARSFT, 
LIFBASE) but many are proprietary. 


| What makes a good measurement technique? S& 


¢ Areliable measurement system vA 
matched with well-understood physics skip 
to make a quantified measurement <a 
that meets a customer’s requirements. ~~ 


Explore each of these... 


| Quantified Measurements 


¢ Accuracy and precision of a measurement 
must be quantified to be useful, for example 
when comparing to a code 
— Qualitative measurements (flow vis) can be useful 


¢ Quantified measurement: 
— Accuracy and precision are well characterized 


Accurate but Precise but ___ (from see ces 
not precise not accurate aeeeny qh Pree len) 


— Accuracy: compared to a ‘standard’ or ‘accepted’ 


— Precision: 10 of large sample of measurements 
¢ Averaging reduces precision’s contribution to uncertainty 
¢ High precision required to quantify fluctuations (turbulence) _ ;, 


| Customer Requirements Ss 


Interview the customer... (sometime the customer is you!) 

¢ Simultaneous measurements to obtain correlations: T’u’, pu’v’ 
¢ Time resolution (< usec), and frequency response (> MHz) 

¢ Spatial resolution: near walls; smallest scales in shear layer 


e Where are the measurement to be made? 
— Inflow to establish B.C., exit, near walls? 
— Point, line, plane, volume? 


e Accuracy and precision requirement? 

¢ What quantity of data is required? Uncertainty required? 

e When is the data needed? Is “real-time” data required? 

¢ What type of optical access is available? 

¢ Can (toxic) seed gases be introduced? Will they influence the 
properties being measured? Can particles be used? 

¢ What is the ordered priority of the above requirements? 


| What makes a good measurement technique? S& 


¢ Areliable measurement system w 4 
matched with well-understood physics (qm 
to make a quantified measurement V4 
that meets a customer's requirements. / 


¢ Now the lectures will focus on different 


measurement systems, introducing the physics. 

— How’s the basic physics of the method work? 

— Can the measurements be made instantaneously 
(aka. single-shot) at high frequencies? 

— What accuracy and precision can be obtained? 


| Customer Requirements S& 


¢ Consider a customer with arbitrary measurement 
requirements: 
— Would like a particular parameter measured 
— Would like the measurement at some conditions (PT) 
¢ Many measurement technique signals scale as Poti Tetatic 
¢ Parameters to be measured: 
— Velocity, temperature, density, pressure, species 
(concentration or mole fraction) 
¢ Can base choice of measurement technique (in part) 
on prior success by others: 


— We compiled past work into graphics that show static 
pressure and temperature of demonstrated measurements 


— Use these charts to find appropriate measurement methods 
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Backup Charts 


| Explanation of Parameter Charts Ss 


Overview: 

This appendix shows the static conditions (temperature and pressure) at which various 
measurement techniques have been successfully demonstrated. When planning a new 
experiment to occur at a specified temperature and pressure, these graphs might be 
useful to see what measurements have previously been demonstrated at those 
conditions. The numerals used as data points indicate the reference in the attached 
reference list. The different measurement techniques are denoted by the colors 
indicated. Pressure is plotted on a logarithmic axis while temperature is on a linear axis. 


Caveats: 


The data includes both references contained in the manuscript as well as additional 
works. 

The reference numbers in the manuscript are different from those in this appendix. 
This study is not meant to be exhaustive; it provides a partial sampling of the data 
available in the literature. 

Some data points overlap. This is particularly true at room temperature and 
atmospheric pressure. Some data has been omitted to avoid overlaps. 

For temperatures greater than 3500 K, the reference number is plotted to the right of 
the chart with the static temperature in parentheses. 

“Species” includes either species concentrations or species mole fraction 
measurements. 


(Not included in AIAA Progress Series manuscript) 


P. M. Danehy, B. F. Bathel, C. Johansen, A. D. Cutler, S. Hurley, “Spectroscopic Measurement Techniques for 


Aerospace Flows,” Von Karman Institute Lecture Series Manuscript, Brussels Belgium, May 12-16, 2014 26 
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Intro to Optical Emission Spectroscopy for ae i ssal 
Hypersonic Nonequilibrium Flows eocien 
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Intro to Optical Emission Spectroscopy for ae i ssal 
Hypersonic Nonequilibrium Flows opine 


Outline: 
¢ Introduction to emission spectroscopy 
¢ Excitation of Atoms and Molecules 
(electronic, vibrational, rotational) 
¢ Line broadening 
¢ Spatial resolution 
¢ Application examples of emission spectroscopy 


Goals: 


te 


Applications of emission spectroscopy ik & 


¢ Astronomy !!! - composition and velocity of celestial bodies 


¢ Diagnostics for gases at high temperatures 
- re-entry (ground test and flight) 
- electric propulsion characterization 
- plasma processing (coating, edging, ...) 
- combustion (though temperatures might be too low) 
- material analysis (laser ablation of surface material > radiation) 
¢ Radiation transport in gases at high temperatures 
- re-entry radiation heat flux to the spacecraft (ground test and flight) 
- radiation transport within the flow field (source term in energy equation) 
> particular importance of VUV radiation 


¢ Identification of elements Typically valid in 
Radiation heat flux determination nonequilibrium 


¢ Temperature determination 
- ratios of atomic lines (electronic excitation temperature) 
- rotational and vibrational temperature of molecules 
(line ratios or comparison with spectral simulation) — 
- translational temperature from line broadening (Doppler effect) Nequilibrium 
- glowing surfaces (Planck blackbody radiation) 
¢ Chemical composition (usually under equilibrium consitions) 


Equilibrium 
assumed but may 
measure thermal 


Typically equilibrium 
required a9 


6) Introduction tk 


> Each particle in a gas has a specific undirected velocity and, therefore, a 
kinetic energy. 

> Definition of translational temperature: for the sum of all particles, a statistic 
distribution function can be formulated (Maxwellian velocity distribution): 


Maxwellian velocity distribution S = 
dN ie V a oe lg 
1 ns ae e “dv “2 


> Similar distribution functions can be formulated 
for other energies such as vibrational, rotational 7 
and electronic excitation (i.e. Boltzmann distributions). 


Energy is constantly exchanged through collisions. 


Electrons are lighter and more movable than heavy particles > more collisons 
possible > possibly different energy distribution > distinctive temperature 

> For each form of energy, a temperature can be defined which 
controls the corresponding energy distribution (Ty, Trots Telecs «++: 


> For many plasmas, the distribution parameters (i.e. temperatures) 
might all be different > T,,ans#Te1 Flo: FT yp DF NON-equilibrium. a 


Emission of Atoms and Molecules zk & 


KENTUCKY 


k=2 A, Vv 
ne OCETTT 
¢ Line strength/intensit hy = 

; Ay, 
¢ Line shape 


eT 

COLT, 
pie 
Co 
PLY i NE 


| 
-10 -8 -6 -4 -2 0 2 4 6 8 10 


from AE=Eypner-Eiower=h Vv; 


> vor Ad for known energy levels 


FWHM 


valid in general (atoms and molecules) 


( hy \_------ number of transitions 
‘ sa nS | > depends on energy of 
A 7 *~~ 1 ~2/ | the upper state 


one photon/transition 
> depends on energy difference 


In equilibrium, the density n, is related to n 


through the Boltzmann distribution: 39 


Emission of Atoms: Electronic Excitation Hs 


KENTUCKY 


¢« Line positio Number of different 


A, Vv 
‘ ition ee | 
: : : states which share the i 
¢ Line strength/intensity same energy 
; Depends on the total 
« Line shape angular momentum J E(A 
which is the value of the 
vector addition of the 
from AE =F upperElower=N a pipial ae ih a 
momentum L an e DNV EW H LN 
> vor a for known energy levels spin momentum S tte ieee 
g=2J+1 


valid in general (atoms and molecules) Degeneracy 


zs opi fede 


hv —s 
Ay) 
Aa - 
Population density of _ 8k E, <— excitation energy of the state k 
; Nn, = VN, €Xp| —— 
the excited state k k U(T. ) 0 T 
ex x ex electronic excitation temperature 
. KL (in equilibrium gas temperature) 
Partition function Total density of the species 


In equilibrium, the density n, is related to n 


under consideration 
through the Boltzmann distribution: 40 


E 
ONES 6 en a 
(7. ) D8 oo 7 


ex 
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¢ Line position 
¢ Line strength/intensity 
¢ Line shape 


from AE=E ner 
> vor Ad for known energy levels 


~Etgwer=h Vv; 


valid in general (atoms and molecules) 


for electronic excitation in equilibrium 


Emission of Atoms: Electronic Excitation x 


CCRC, 


eee lee 
fede 


information on the transition (from 
quantum mechanics, e.g. from tables) 


information about the thermodynamic 
condition of the plasma 


In equilibrium, the density n, is related to n 
through the Boltzmann distribution: at 


tk & 


Li ' . . 
inf +const = — Ey > Plot nf to vs E, > straight line with slope 1/-kT 
VA, 2; kT... (Boltzmann Plot) 


Electronic Excitation Temperature 


Theoretically, the ratio of two lines would be sufficient to determine the excitation 
temperature, but it is often in doubt if the electronic excitation is in equilibrium. 


If points are indeed on a straight line > Boltzmann valid > T can be determined 


BUT: Even if Boltzmann seems valid, the apparent excitation temperature might not be an 
equilibrium temperature: 


Fa) “5 ——S 
{e) 

1S) 

+ -10 

“Zz 

ed 

= 

= -15 ~S3=76000K 

£ ena decreasing T 

© -20 


=== 128000 K 


“== 256000 K 


-25 
95000 97500 100000 102500 105000 107500 110000 


electronic excitation energy E,, cm 
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Ws 


Electronic Excitation Temperature tk Gy 


Theoretically, the ratio of two lines would be sufficient to determine the excitation 
temperature, but it is often in doubt if the electronic excitation is in equilibrium. 


Li ' . ' 
inf +const = — Ey > Plot nf to vs E, > straight line with slope 1/-kT 
VA, 2, (Boltzmann Plot) 


If points are indeed on a straight line > Boltzmann valid > T can be determined 


ex 


BUT: Even if Boltzmann seems valid, the apparent excitation temperature might not be an 
equilibrium temperature: In(n/g) 


recombining plasma, 
_ electronic overpopulation 
“hy UXTP UPL. 1 lee 


Recombining plasma 
— high energy states will be populated through 
recombining ions 

> flat distribution 

> high excitation temperature 


lonizing plasma 

— high energy states will be depopulated 
trough ionizing neutrals 

> steep distribution 

> low excitation temperature 


ionizing plasma, 
electronic underpopulation 
LTA a ToT. 


Emission from Molecules Ws 
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Molecule emission shows 


some shape but looks like a 

Sos mm StonrerssooK continuum 

= Trot=8100K = 

5 06 - Ratio of the different vibrational 

Bod systems governed by T,,, 

ee 

3 0.2 ¢ Ratio within one vibrational transition 

(e.g. 0-0 > V,=0, V|=0) governed by T,; 

0 


300 320 340 360 380 400 420 
wavelength, nm 


0-0 
“a om sim. 
High spectral resolution es a 
>Single lines like atoms " 1-H = 8100 Monelelsiieaice 
ut more ... ". 
0,3 


> rotation and vibration 


relative intensity |-| 


) Ae UAV 


388 388,5 389 389.5 390 390.5 39] 3915 
wavelength, nm 
Winter, M.: Emissionsspektroskopische Untersuchung der Umstr6mung von Probenkérpern in hochenthalpen Plasmastr6mungen / Emission Spectroscopic Investigation of the 
Flow Field around a Blunt Body in a High Enthalpy Flow, Dissertation (in German), Universitat Stuttgart, Germany, Dec. 2006. 44 


Emission from Molecules ik 


Radiation of molecules is caused by changes of 
the electronic level, the vibrational level, and the 
rotational level of excitation, often 
simultaneously. 


Born-Oppenheimer approximation: The total energy is the sum of the 
different excitations which can be computed separately. 


Electronic excitation works similar to the atoms ... 
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Emission from Molecules es 


KENTUCKY 


¢ For diatomic molecules, rotation and only one vibrational mode need to be accounted for. 
¢ Already for tri-atomic molecules, several vibrational modes are possible. 


I : ; 
diatomic 
—§* molecules 
I 
linear triatomic linear triatomic asymmetric « Potstionandwibration 
symmetric stretch mode stretch mode 


energies are discrete and not 


continuous and are obtained 
-— from solutions of the 


Schroedinger equation. 


e—-e-6 


¢ Different modes for linear and non-linear molecules 


¢ For diatomic molecules, each 
energy is characterized 
through one quantum number 
(J for rotation, V for vibration). 


30 : "2 
symmetric stretch asymmetric stretch 


a V4 \ 


¢ In addition, molecules may be 
b> electronically excited 


. (Same process as for atoms). 
vibration of non-linear triatomic molecules 


Main vibrations occurring in liquid water 
(http://www .lsbu.ac.uk/water/water_vibrational_spectrum.html) 46 


Emission from Molecules 


* Each electronic level of a molecule is described by a potential 
curve in terms of energy vs. inter nuclear distance. 


upper electronic 


rotational energy level E.. or T 
el,2 2 


ELota(Jo) OF Froto(Jo) 


vibrational energy 
Evin,a(Va)Or Gyipo(Vo) 


lower electronic 
level E,,, or T, 


internal ener 


rotational energy 
Eot(d1) or Frot(dy 


\..... 2=/ V:=2 vibrational energy 
a Evin(V1) Or Gyin(V4) 


intranuclear distance 


Vv = + (AE, + AE, ip + AE,,,) 
he 


tk & 


* Different vibrational 
energies appear as lines of 
constant energy ion these 
potentials 


* Excitation of the levels can 
happen through collisions 
(one major partner are 
electrons) or absorption of 
radiation. 


* From one upper level, 
transitions to different lower 
levels are possible, each 
one resulting in a different 
emission line 


¢ Usually, rotational energy 
changes in the same 
transition, possibly electronic 
excitation does, too, 
therefore contribution to the 
transition energy and the 
line position. 


47 


Emission from Molecules 
¢ Each electronic level of a molecule is described by a potential 
curve in terms of energy vs. inter nuclear distance. 


upper electronic 


level E.,, or T, 


V9o=—G4 


lower electronic 
level E,,, or T, 


Energy 


inter nuclear distance 9——~> 


Vibrational Excitation of diatomic molecules 


tk & 


* Different vibrational 
energies appear as lines of 
constant energy ion these 
potentials 


¢ Excitation of the levels can 
happen through collisions 
(One major partner are 
electrons) or absorption of 
radiation. 


* From one upper level, 
transitions to different lower 
levels are possible, each 
one resulting in a different 
emission line 


¢ Usually, rotational energy 
changes in the same 
transition, possibly electronic 
excitation does, too, 
therefore contribution to the 
transition energy and the 
line position. 
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Emission from Molecules 
¢ Each electronic level of a molecule is described by a potential 
curve in terms of energy vs. inter nuclear distance. 


upper electronic 


level E.,, or T, 


= - 7 VY 


Vo=—G 


lower electronic 
level E,,, or T, 


Energy -—— 


———==> 


inter nuclear distance 


Vibrational Excitation of diatomic molecules 


tk & 


* Different vibrational 
energies appear as lines of 
constant energy ion these 
potentials 


¢ Excitation of the levels can 
happen through collisions 
(One major partner are 
electrons) or absorption of 
radiation. 


* From one upper level, 
transitions to different lower 
levels are possible, each 
one resulting in a different 
emission line 


¢ Usually, rotational energy 
changes in the same 
transition, possibly electronic 
excitation does, too, 
therefore contribution to the 
transition energy and the 
line position. 
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Emission from Molecules tk 


Rigid Rotator and Harmonic Oscillator 

- simple model, mainly to be used to account for 
rotational and vibrational energies in CFD codes > 
energy conservation 

- not suitable (at least the harmonic oscillator for 
spectral simulation) 

¢ The harmonic oscillator, however, does 

not describe real molecules: 

¢ forr~inf. > E =dissociation energy 
*forr~O D> E~ int. 


\ q i harmonic 
nat iS i oscillator 
— Y ! 

o is / 

Cc = / 

® = / 

© a | dissociation energy 

¢ = 

Sam joy 

® q 

= - . . 

= l realistic 
potential 


intra-nuclear distance 


Morse Potential 

U(x) =D,(l-e* J 

* x: distance to r, 

(re: equilibrium distance between 
the nuclei) 


* D,: Well depth of the potential 
(dissociation energy- energy at 
minimum) 


B =1,2177-107 @, |“ 
D 


e 


* ua: reduced mass of the molecule 


vibrational energy increment from V 
to V+1 now decreases with V: 


G(V )=a, [v + 7 — 0X, iv + 7 
Z 2 
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Emission from Molecules es 


KENTUCKY 


¢ Separation of rotational and vibrational excitation is not completely possible. 
¢ In fact, the whole molecular potential changes shape with rotational quantum number. 
* For rotationally highly excited states, only low vibrational excitation possible > V_,,,(J) 
¢ Or: limiting rotational quantum numbers different for each vibrational state > J,.,(V) 
Molecular potential with rotational excitation 
200000 ne ©) * limiting case: 
= no more potential well 
> not stable 


¢ will be different for each 
electronic state 


180000 
160000 
140000 i a 
120000 
100000 
80000 


energy, cm"! 


60000 
40000 
20000 


0 
0,05 01 0,15 0,2 0,25 0,3 0,35 
inter nuclear distance, nm 


Winter, M.: Emissionsspektroskopische Untersuchung der Umstrémung von Probenk6rpern in hochenthalpen Plasmastr6mungen / Emission Spectroscopic 
Investigation of the Flow Field around a Blunt Body in a High Enthalpy Flow, Dissertation (in German), Universitat Stuttgart, Germany, Dec. 2006. 51 


6) Emission from Molecules tk 


Theoretical Simulation of Molecule Radiation 


> Intensity of one emission line: 


_ NA, AE., ~ Ly \ l6m°cv" 2 ae 
a ioe: ee or Se VN} ees R _ _ J"A 
An ~ 3 enn av 2s'+1 
—~ vW 
Particle density Electronic Honl London Factor 


in the levele’ V’' J’ transition moment 


esunipiion et N 0%0"" of! with the total aes function: 
ev! dibs ror _ ov 
Boltzmann distributions = Niji ==>, Of Of Oey? o- » 20.0; 
evs’ . 
a g, (2/'+1) E‘ Ey Eo 
N ‘j i tr, IN ——. 7 © : ex = exp| — —?-— 
q oe Ta J Nm QO (er) kt) KD, 
electronic vibrational rotational 


e' ES 'V' Eo Vis! E 
= &,, exp| ——— = us : 2s'+1 a ae 
Q., § e w kT Qin = exp [- Over = ( Jexp kT 


el k Dip 
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Emission from Molecules Ws 
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Theoretical Simulation of Molecule Radiation 


o= =~ 


7 N 
> Intensity of one emission line: rs 
' 3-4 / \ JOA 
ae N ou _WN l6z°cv (R, Foye) dye San 
1 > « 2J'+1 
| \————"’ 1 Da 
Particle density‘. Electronic ,/ Hodnl London Factor 
in the level e’ V’ J’ transition moment 

* The electronic transition is fast in comparison Frank-Condon Principle 


to the motion of the nuclei during vibration. 


>The position of the nuclei relative to each 
other will not change during the transition. 


b) 


Condon factor. 


Figure from Demtroder, Wolfgang. "Atoms, molecules and photons." Heidelberg, 
Springer-Verlag Berlin Heidelberg, 2006. 571 (2006) with permission. 


> The higher the overlap between the upper Epot(R) 
and lower potential, the higher the probability of 
the transition. 
> If wave functions are analyzed, this ies 
probability can be expressed in the electronic 
transition moment containing the Franck- ! os 


Vs 


KENTUCKY 


Molecular Band Structure 


¢ Transition energies with the same difference in vibrational quantum number Av are similar 
> The different vibrational transitions tend to group for the same Av. 


*AV =0 : 0-0, 1-1, 2-2, 3-3, ... 


sau a upper level e’ 
NOy 
Tiot=2500K p 0,045 
— 2,5E-05 5 Tyib=10000K 
7 : 0,04 > 
7) @ oD) 
£ 7" _ easurement L 0,035 Ss ® 
© 2,0E-05 - — EMCAL simulation . Cc 
= 2 Ww 
£ 0,03 
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Michael W. Winter and Monika Auweter-Kurtz: Status on Molecular Simulation Capabilities of PARADE and EMCAL and Comparison with Experimental 
Data in High Spectral Resolution, 5° European Symposium on Aerothermodynamics for Space Vehicles, KéIn, Deutschland, 8—11.11.2004. 54 
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Molecular Band Structure 


¢ Transition energies with the same difference in vibrational quantum number Av are similar 
> The different vibrational transitions tend to group for the same Av. 


*AV=1: 1-0, 2-1, 3-2, 4-3, ... > higher AE > lower A 


ais a upper level e’ 
NOy 
Trot=2500K p 0,045 
2,5E-05 + Tyib=10000K 
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ey 
— measurement — |L 0.035 © 
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2S 
oO 
ine) 
o1 
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Michael W. Winter and Monika Auweter-Kurtz: Status on Molecular Simulation Capabilities of PARADE and EMCAL and Comparison with Experimental 
Data in High Spectral Resolution, 5° European Symposium on Aerothermodynamics for Space Vehicles, KéIn, Deutschland, 8—11.11.2004. 95 
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Molecular Band Structure 


¢ Transition energies with the same difference in vibrational quantum number Av are similar 
> The different vibrational transitions tend to group for the same Av. 


© AV =-1 : 0-1, 1-2, 2-3, 3-4, ... > lower AE > higher A 
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Michael W. Winter and Monika Auweter-Kurtz: Status on Molecular Simulation Capabilities of PARADE and EMCAL and Comparison with Experimental 
Data in High Spectral Resolution, 5° European Symposium on Aerothermodynamics for Space Vehicles, KéIn, Deutschland, 8—11.11.2004. 56 
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Temperature determination from molecular radiation 


¢ Within one vibrational transition, a large number of rotational lines is present 
> rovibrational lines (may shift to lower or higher wavelengths, depending on B) 


* Similar to the atoms, a Boltzmann plot can be performed to determine T,,;. 


wavelength, nm 
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Probe in a Subsonic Plasma Flow,” 31st AIAA Thermophysics Conference, New Orleans, LA, June 17-20, 1996. 


| Abel-Inversion KENTUCKY 


— emission spectroscopy is a line of sight method 
: plasma 
> only integrated data, no spatial resolution ealinne 


focusing 
lens 


\ detector 


Abel integral equation: 


I(z) = {SS 


e(r) dr 


Typically, integrals along the line of sight are measured 
> the local emission is desired. (i 
ee 


whe Jon? 


dz 
Solving this equation is usually only possible in 12 —= yal 
special cases (e.g. a constant local distribution 
e(r)=const. yields and elliptical profile in I(z)) 
Under the assumptions of rotational symmetry 


and optically thin medium, local emission values —; ey 
can be obtained from an Abel-Inversion. 
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Approximate Abel-Inversion ek 


2.5E+09 8.0E+07 


> Different radial/vertical 
positions are to be 
measured 


> Amatrix can be built 
that gives a relation 
between integrated 
and local emission 


2.0E+09 = 6.4E+07 


1.5E+09 4.8E+07 


1.0E+09 rw 3.2E+07 
—*— NEQAIR, int vu 13 


5.0E+08 -—<—vu}t+3tocat 1.6E+07 


local population density, cm? 


LOS-integrated population density, cm 


0.0E+00 = t t 0.0E+00 
0 50 100 150 200 


distance off-axis, mm 


> The matrix simplifies significantly ifonlyfew = ------- 
measurement positions are taken into account. 
constant local intensity in each ring assumed 

> Avery simple system of linear equations can be obtained 
and solved recursively 


Is int (Z3) = e(r;) *L,(r3) 

lint (Z2) = e(rz) *L,(r3)+ — e(r,) *L,(r,) 

I int (21) = e(r3) *L,(r3)+  e(r,) *L,(r,)+ e(r,) *L,(r,) ZZ Z, Z, 
lo,int (Zo) = e(rs) *Lo(r3)+ el.) *Lo(r,)+ elr,)*L(r)+ el) 


Winter, M. W., Prabhu, D. K., "Excited State Chemistry in the Free Stream of the NASA IHF Arc Jet Facility Observed by Emission 
Spectroscopy," 42nd AIAA Thermophysics Conference in Honolulu, Hawaii, 27 - 30 Jun 2011. 59 


Approximate Abel-Inversion 


2.5E+09 8.0E+07 


2.0E+09 = 6.4E+07 


positions are to be 
measured 


> Amatrix can be built 
that gives a relation 
between integrated 
and local emission 


1.5E+09 4.8E+07 


1.0E+09 rw 3.2E+07 
—*— NEQAIR, int vu 13 


5.0E+08 -—<—vu}t+3tocat 1.6E+07 


local population density, cm? 


LOS-integrated population density, cm 


0.0E+00 = t t 0.0E+00 
0 50 100 150 200 


distance off-axis, mm 


> The matrix simplifies significantly if only few 
measurement positions are taken into account. 
constant local intensity in each ring assumed 


> Avery simple system of linear equations can be obtained 


and solved recursively 


Is int (23) = e(rz) *L,(r3) 

lp int (Zp) = e(r3) *Lj(r3)+ —e(r,) *L,(r,) 

ly int (24) = e(r;) *L,(rz)+  e(r,) *L,(r,) +  e(r,) *L,(r,) 
lo int (Zo) = e(r;) *Lo(r,)+  e(r,) *L(r,)+  e(r,) *L,(r,) + 


> Different radial/vertical 
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Winter, M. W., Prabhu, D. K., "Excited State Chemistry in the Free Stream of the NASA IHF Arc Jet Facility Observed by Emission 


Spectroscopy," 42nd AIAA Thermophysics Conference in Honolulu, Hawaii, 27 - 30 Jun 2011. 
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Application examples chien 


Measurements in an atmospheric, inductively coupled plasma torch 
( C. Laux, Stanford University, now Ecole Centrale, Paris) 


¢ Equilibrium plasma at up to 6000 K 
> chemistry well defined 


¢ Used for the development of radiation models * After equilibrium characterization, 


used for non-equilibrium investigation. 
¢ Excellent agreement between experiment 


and simulation can be achieved. 


a or LO 7, 0.20 
| —— SPECAIR: NO B,,3,¢,8',7' + O, SR + atomic lines ™ i N, (B>X) 
Q (with self and room air absorption) Pos 08 LL - Measured 
Mg a 7 § t N. (CB) —— SPECAIR-LTE 0.15 
oO ao L + 
2 aes 
Boat & r 0.10 
8 i+ 0.4 i N, (B —> A) 

& Z 2 ; 
a *T, 2 02F 0.05 
oL ll 0 ~ 7 . ~ - 0 

300 350 400 450 600 650 700 750 800 


: TT Pe ee ee ey ee ae ae 

140 160 180 200 220 240 260 280 
Wavelength (nm) Wavelength (mm) 

Laux C. O.: Radiation and Nonequilibrium Collisional-Radiative Models, VKI Special Course on Physico-Chemical Models for High Enthalpy and Plasma 
Flows Modeling, Belgium, June 4-7, 2002. See also Laux, C.O., Pierrot, L., and Gessman, R.J., “State-to-state modeling of a recombining nitrogen plasma 
experiment,” Chemical Physics, Vol. 398, pp. 46-55, DOI: 10.1016/j.chemphys.201 1.10.028, 2012. 
Laux, C. O., “Radiation and nonequilibrium collisional-radiative models” in “Physico-chemical Models for High Enthalpy and Plasma Flows, VKI LS 2002-07, 
edited by D. G. Fletcher et al., Rhode-Saint-Genése, Belgium, 2002. 61 
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Winter, M, Srinivasan, C., Charnigo, R., “Non-Equilibrium Analysis of Emission Spectroscopy Data Taken in the Freestream of ine 
NASA IHF Arc Jet Facility,” AIAA AVIATION 2015, Dallas, Texas, June 2015. 
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Comparison of free stream 
spectra at high condition with 
LTE simulation (VIS/NIR): 


¢ Spectra dominated by 
emission of N2 B-A (1st Pos.) 

¢ Thermodynamic quantities 
from CFD solution used as 
input for NEQAIR 

¢ Spectra integrated along 
the line of sight 


¢ Simulation and experiment 
do not even agree qualitatively 


¢ If individual vibrational states are 
computed separately, the spectrum 
can be fitted by scaling these upper 
populations 

> massive overpopulation of the 
uo vibrational levels peaking at 


Vupper = 13 
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Inverse level crossing with N, BSI1, at v between 10 and 13 
Predissociation 


> Population of the excited state is not 
caused by thermal excitation of the 
ground state but by recombination of 
atoms in an excited N, state 


cen, 


bags 
2 
~ 1 . 
- > Molecule emission can be used to 
: determine atom ground state densities. 
g <pen 
2 2E+16 
ui CFD 
E eons @| vu=12, Benson rate, Acton 
a 1-6E+16 > vu=13, Laux rate, Acton 
2 1.4£+16 = =13, ; 
¢ 1.2E+16 -::-:= --- 8A 0 -vu=12, Benson rate, Ocean Optics — 
mo} E416 eee = - 
a. 7 SE+15 Lee 
i ~xtyt & wet 
: sd © GE+5 amen 
ots SRNEK [ee ee Ee! OOO ee oe % wo 
2.5 3.6 4.5 55 ° 4E+15 TS, 
BOND LENGTH, a, = 2E+15 
¢ Observed radiation comes from the 1! Pos. 0 
System (N, BI, > A %z,* ) with electronic : : - = 20 
‘ i 9g radial distance, cm 
excitation energies around 8eV to 10eV. 
Harry Partridge, Stephen R. Langhoff, and Charles W. Bauschlicher, Jr., and Winter, M, Srinivasan, C., Charnigo, R., “Non-Equilibrium Analysis 
David W. Schwenke, “Theoretical study of the A°S,* and C’”®P,, states of No: of Emission Spectroscopy Data Taken in the Freestream of the 
Implications for the N, afterglow,” J. Chern. Phys. 88 (5), 1 March 1988. NASA IHF Arc Jet Facility,” AIAA AVIATION 2015, Dallas, Texas, 
Figure reproduced with permission from the authors. June 2015. 63 


Application examples 


Position (cm) 


Mean Spectral Radiance 
(W/cm3-sr-um) 


Quantify radiation loads to spacecraft during 
atmospheric entry in impulse facilities 


Shot 35 (Blue) - 8.53 km/s, 0.098 torr 
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192.20 to 345.85 nm 


Goals: 
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Shock 
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EAST shock tube 
(NASA Ames): 


* reproduction of flight 
conditions 

* no model testing 

¢ 130nm — 8um in various 
intervals 

* simultaneous measure- 
ments with 4 spectrometers 
possible 

* one dimension resolves 
spectrally, the other spatially 
> shock and post-shock 
region are covered in one 
image 


* quantify radiation from shock and post shock 

> radiation heat flux from the plasma upstream of 
the spacecraft after spectral integration 

> empirical correlation for radiation heat flux 


- Cruden, B., Martinez, R., Grinstead, J. H., Olejniczak, J., "Simultaneous Vacuum Ultraviolet through Near IR Absolute Radiation Measurement with 
Spatiotemporal Resolution in an Electric Arc Shock Tube," AIAA 2009-4240, 41st AIAA Thermophysics Conference, 22 - 25 June 2009, San Antonio, Texas. 

- Aaron Brandis, Christopher Johnston, Marco Panesi, Brett Cruden, Dinesh Prabhu, Deepak Bose, “Investigation of Nonequilibrium Radiation for Mars Entry,” 
AIAA 2013-1055, 51st AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, 2013. 64 
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Application examples 


Monitoring of ablation species through optical methods 


Testing of PICA at NASA Ames (SPRITE) showed the appearance of 
ablation products in emission spectra inside the boundary layer: 
- CN, OH, NH, Ca, (Na, kK) 


Combination of pyrolysis products themselves 
and interactions with the plasma. 


[oe] 


W m2? sr? nm 


spectral radiance, 


300 320 340 360 380 400 420 440 460 480 500 
wavelength, nm 


D. Empey, K. Skokova, P. Agrawal, G. Swanson, D. Prabhu, K. Peterson, M. Winter, E. Venkatapathy, “Small Probe Reentry Investigation for 
TPS Engineering (SPRITE),” AIAA-2012-0215, 50th AIAA Aerospace Sciences Meeting, Nashville, Tennessee, 9 - 12 Jan 2012.. 65 


Application examples zk & 


Remote recession measurements of ablating TPS material 
LJ Tracer element (Coating/paint) (on ground and in flight) 
Tracer plug (e.g. PICA) 
CN 


l F 


Less recession than 
seeding depth 


Michael W. Winter, Margaret Stackpoole, Anuscheh Nawaz, Gregory Lewis Gonzales, and Thanh Ho, “Remote Recession Sensing of Ablative 
Heat Shield Materials,” 52nd AIAA Aerospace Sciences Meeting, 13 - 17 January 2014, National Harbor, Maryland. 66 


Recession exceeds 
seeding depth 


Recession reaches 
seeding depth 


Conclusions x & 


Optical emission spectroscopy is a useful tool to access information on thermal 
characteristics and composition of a plasma. 


Applications for which the assumption of equiliorium is rather uncritical: 
- identification of species; 
- total radiation flux (in the wavelength covered and if properly calibrated); 
- quantities from line broadening (e.g. Doppler temperature, electron density); 
- identification of thermal non-equilibrium. 


Applications for which the assumption of equilibrium is rather crucial: 
- determination of particle densities from thermal excitation; 
- measurement of excitation temperatures from line ratios. 


Proper calibration is imperative. 
Comparison with simulation (non-equilibrium chemistry and excitation, e.g. 


though collisional radiative models) seems the best approach as soon as these 
models are validated. 


Thanks for your attention! Open for 


Questions ? 
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Back-up Slides xk 


Selection of Simulation Codes (Line-byLine) 


"Fx 


NEQAIR (Non Equilibrium ...) - NASA Ames EAR restricted 

FORTRAN, non-eq. (QSS + multi T), air and CO, plasma related species + H, radiation transport capabilities 
Whiting, E E, Park, C, Liu, Y, Arnold, J O, and Paterson, J A, “NEQAIR96, Nonequilibrium and Equilibrium Radiative Transport and 
Spectra Program: User’s Manual,” NASA RP-1389, NASA, December 1996. 

HARA (...) - NASA Langley EAR restricted 

FORTRAN, non-eq. (QSS + multi T), air and CO, plasma related species, radiation transport capabilities 
Johnston, C. O., Hollis, B. R., and Sutton, K., “Spectrum Modeling for Air Shock-Layer Radiation at Lunar-Return Conditions," 
Journal of Spacecraft and Rockets, Sep.-Oct. 2008. 

SPRADIAN/RADIPAC — Japan/South Korea available on request to JAXA 


FORTRAN, non-eq. (QSS + multi T), air + ablation related species , radiation transport capabilities 
Fujita, K., and Abe, T., “SPRADIAN, Structured Package for Radiation Analysis: Theory and Application,” 


The Institute of Soace and Astronautical Science Report No. 669, July 1997. 
PARADE — European Space Agency ESA available on request to ESA 
FORTRAN, non-eq. (QSS + multi T), air + CO, plasma related species , rad. transport capabilities (HERTA) 


Smith, A. J., Wood, A., Dubois, J., Fertig, M., Pfeiffer, B.: Plasma Radiation Database PARADE V22 Final Report Issue 3, 
ESTEC contract 11148/94/NL/FG, FGE TR28/96, Issue 3, Oktober 2006. 


SPECAIR — Ecole Centrale Paris Windows version available by download 
air, multi T, no radiation transport 


Laux, C. O. , "Radiation and nonequilibrium collisional-radiative models” in “Physico-chemical Models for High Enthalpy and 
Plasma Flows,” VKI LS 2002-07, edited by D. G. Fletcher et al., Rnode-Saint-Genése, Belgium, 2002. 


http://www. specair-radiation.net 


LIFBASE - Stanford Research Institute SRI Windows version available by download 


LIF related diatomic molecules (CH, NH, OH, CO, NO, ...), no radiation transport 
J. Luque and D.R. Crosley, "LIFBASE: Database and spectral simblation (version 1.5)", SRI International Report MP 99-009 aed, 
http://www.sri.com/psd/lifoase/ 


Intensity Calibration — How? ek & 


5000 


Integrating Sphere 


Deuterium Lamp Calibration Lamps 


— Deuterium Lamp 


Deuterium lamp / uW/(cm? sr nm) 
Integrating Sphere / wW/(cm? sr nm) 


1000 — Integrating Sphere 50 
0 0 
200 400 600 800 1000 


wavelength / nm 


¢ Calibration with an integrating sphere ( halogen bulbs - UV-VIS-NIR) and a Deuterium an (UV) 
both calibrated to spectral radiance in uW/(cm? sr nm) 


* Correction factor: ratio of factory calibration value and measurement of calibration lamp 


¢ Deuterium lamp (UV) can not easily be used for absolute calibration 
> scaling of Deuterium correction factor to the |-Sphere in overlapping region. 

¢ Best calibration: place the calibration lamp at the location of the measured plasma 
> all set-up influences are covered and direct calibration is possible. 


¢ Lamp area must be larger than measurement spot !!! 


Winter, M. W., Prabhu, D. K., "Excited State Chemistry in the Free Stream of the NASA IHF Arc Jet Facility Observed by Emission 
Spectroscopy," 42nd AIAA Thermophysics Conference in Honolulu, Hawaii, 27 - 30 Jun 2011. 70 


Intensity Calibration — How? 
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ee £OPFneg Eoew meas (4) 


ts 


KENTUCKY 
10 
get 
” 
G18 — Corr. Sphere 
c 
a 1-6 ——Com Deut 
-- | < 1.4 (Scaled) 
E 1.2 
G 4 
5 0.8 
- 01 306 
_ 0.4 
Fr 0.2 
° 
O o 
0.01 320 340 360 380 400 
4000 wavelength / nm 
_ dp tisees (A) I A 
~~ I (A en | ) 
BB,meas ) 
dt ene (A) 


> Correction factor corr(A) which is the inverse sensitivity of the set-up. 


> Use of Deuterium lamp only for qualitative calibration 


followed by scaling (cross calibration) through halogen lamp/black-body. 


Winter, M. W., Prabhu, D. K., "Excited State Chemistry in the Free Stream of the NASA IHF Arc Jet Facility Observed by Emission 
Spectroscopy," 42nd AIAA Thermophysics Conference in Honolulu, Hawaii, 27 - 30 Jun 2011. 7 
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Emission of Atoms: Electronic Excitation ik & 


hv hv g E 
— AG A; nm = an ory : 


4a 
¢ For low temperatures, only the first energy levels are populated. 
¢ With temperature going to infinity, all levels would be populated according to g,/U. 
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Electronic Excitation of Hydrogen tk 


E 
A, n, = —A, n,exp| ——<* 
ok Ax 1 Gee a 


normalized line emission €,, 


450 500 550 600 650 700 
wavelength, nm 


Im?,Wm? srt 


emission &,, for ny 


Electronic Excitation of Hydrogen zk & 
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Balmer spectrum of hydrogen for different temperatures: 
(logarithmic scale) er 

¢ Emission increases by 60 orders of 
magnitude between 1000K and 8000K 


—— 128000 K 
—— 64000 K 
—— 32000 K 


¢ Ratio of different lines in the Balmer 
series changes with temperature 
(hard to see on log scale). 


—— 16000 K 


—— 8000 K 

—— 4000 K ¢ To illustrate the ratios, we normalize 
—2000k | the spectrum to the H, line. 

= 1000 K 


wavelength, nm 
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Electronic Excitation of Hydrogen 


Ki Ny 


normalized line emission &,, - 


tk & 


Spectrum normalized to Ha to illustrate 


the influence on the line ratios: 


300 350 400 450 500 550 600 650 700 
wavelength, nm 
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Emission of Atoms and Molecules tk 


Line Broadening 


Emission and absorption lines can be broadened through different processes: 

¢ natural broadening (Heisenberg uncertainty principle) 

¢ Doppler broadening (thermal motion) 

* Collision broadening 
- Van-der-Waals: collisions with neutrals 
- Resonance: collisions with the ‘like’ species particles (perturber’s state connected by 

an allowed transition to the upper or lower state of the transition under consideration) 

- Stark for interactions with ions or electrons, specific for each transition/line 

¢ (Zeeman effect - interaction with magnetic fields) 


¢ Each of these processes is actually 
a line shift for one photon. 
¢ The integration over many photons with 
different shifts produces a broadening. 
> distributing the transition energy over a 
certain wavelength/wavenumber range 


¢ Depending on the physical effect, the final line shape will 
follow a Gauss or a Lorentz profile. 
¢ Usually, the broadening is described by the line width at 


i © Half widths of Lorentz profile can be summed up linearly. 


Wr 
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Line Broadening 


Collision broadening (Lorentz) Doppler broadening (Gauss) 
S be f 
K;, = ee’ S = K; dn, — 
"2am? +B an = 


HWHM 7 (2) 
decreases with T Po 


Lorentz: 

* 90% of the transition 

energy within 6 HWHM 

* 99% is reached after 
~ 50 HWHM 


Gauss: 
° 99% in 2 HWHM 


fraction of transition energy 
inbetween2xAn 


An/HWHM, - 
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The Anharmonic Oscillator — Higher Order Potentials 


Energy [cm] 


* Higher order potentials 


ie Morse Potential 
Lippincott — 3 parameter fit 


Hulbert-Hirschfelder — 5 parameter fit U(x) = D, ( _ eh 
¢ may provide better agreement with 
experimental data * x: distance to r, 
(re: equilibrium distance between 
1.E+05 the nuclei) 


Hulbert-Hirschfelder and 


9.E+04 Morse Potentials 


8.E+04 
7.E+04 
6.E+04 
5.E+04 
4.E+04 
3.E+04 
2.E+04 
1.E+04 
0.E+00 


* D,: Well depth of the potential 
(dissociation energy- energy at 
minimum) 

B=1,2177-10" a 


e 


* u,: reduced mass of the molecule 


NOB H.H. NOC H.H. 
NOD HH. g-=enox mu. Vibrational energy increment from V 


—=—=NOAM- NOB M. to V+1 now decreases with V: 
NOC M. NOD M. 


1 1) 
0.5 1 1.5 2 2.5 3 av)=olv +5 |-or{v+s) 


inter-nuclear distance [A] 
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Higher Order Approximations - Dunham Expansion 


¢ Vibrational energy expressed as a polynomial function of V 
¢ For small quantum numbers equivalent to Morse 


l 1) 1) 1) 
GV )=0, V+— |-@,x,; V+—| +o,y,|; V+—]| +@,z,| Vt—] +... 
2 2 2 2 


¢ Rotational energy contains higher orders of J 
* coupling constants now depend on V 


F,(V,J)=B,J(J +1)—D, J? (J +1) +... 


B.=B-a vel + vel 7 

V e e 5 Va 9) eee 
1 1yY 

D, =D,+£. es PO\ Vl 


¢ Different sets of constants given in literature 


Dunham, J: The energy levels of a rotating vibrator, Phys. Rev. 41, 721, 1932. 
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¢ Best calibration: place the calibration lamp at the location of the measured plasma 
> all set-up influences are covered and direct calibration is possible. 


¢ Lamp area must be larger than measurement spot !!! 
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Intensity Calibration — How? lk & 


1200 ery 10 


UV-VIS-NIR = =Deut measured 
— Deut., calibrated 
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(jo) 


Corr., combined 


Oo 
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radiance / pW/(cm? sr nm) 
measured intensity / counts/ms 
(op) 
(je) 
oO 
Corr. / pW/(cm? sr nm)/(counts/s) 


400 N s- 0.1 
200 l 
I 
0 oe — ia 0.01 
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wavelength / nm 
¢ Deuterium lamp (UV) uses a hydrogen discharge to produce high radiation in VUV-VIS 


¢ Usually, only the continuum emission from UV down is used for calibration. 
¢ Although calibrated to radiance in uW/(cm?2 sr nm), the dimensions of the discharge are small 
(~1mm) > measurement spot usually larger 


¢ Discharge shows gradients > additional uncertainties introduced. 


> Use of Deuterium lamp only for qualitative calibration 
followed by scaling (cross calibration) through halogen lamp/black-body. 


re 


Wavelength Calibration tk & 


¢ Compare measured line positions (pixel) with known line positions (e.g. Hg lamp): 
- scanning spectrometer: typically linear with time > at least 2 lines needed; 
- CCD measurement: typically binomial > at least 3 lines needed. 

¢ additional lines may serve as a calibration check ... 
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Electronic Excitation Temperature "rx 


KENTUCKY 
The emission coefficient is integrated along the line of sight 
hv 
pa A n Te =| eQ)dx=—— Aynyl 
Arr 47 
With the Boltzmann distribution follows: 
E I. E 
n= Ok No exp - a in a +const = ——4 
U(T.,) kT. VA 8 KT. 
0 


Theoretically, the ratio of two lines would be 
sufficient, but it is often in doubt if the 


an eS 
electronic excitation is in equilibrium. c 
>) 
-10 
> Boltzmann Plot: ae 
on 
I, = 
Plot nf re VS E,. ye 15 += 76000K 
ee “3 =—@—32000K decrea sing T 
> Straight line with slope 1/-kT = 20 Ste Pel 
=== 128000 K 
If points are indeed on a straight line -25 (286000 K SS 


+> Boltzmann valid 95000 97500 100000 102500 105000 107500 110000 
+> T can be determined electronic excitation energy E,, cm™* 
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Molecular Band Structure Ue 
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¢ Transition energies with the same difference in vibrational quantum number Av are similar 
> The different vibrational transitions tend to group for the same Av. 


© AV =-1 : 0-1, 1-2, 2-3, 3-4, ... > lower AE > higher A 
3,0E-05 0,05 upper level e’ 
i y 0,045 
— 2,5E-05 5 Tyip=10000K 
= L 0,04 


2,0E-05 + —EMCAL simulation 
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— measurement 
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ine) 
o1 
simulated emision [arb. units] 


1,0E-05 + 
r 0,015 
0,01 
5,0E-06 
Av=2, 
| Aves : oe 0,005 
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vacuum wavelength [nm] 


AV=-1: transition energy lower than AV=0 


Michael W. Winter and Monika Auweter-Kurtz: Status on Molecular Simulation Capabilities of PARADE and EMCAL and Comparison with Experimental 
Data in High Spectral Resolution, 5° European Symposium on Aerothermodynamics for Space Vehicles, KéIn, Deutschland, 8—11.11.2004. 84 
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Theoretical Simulation of Molecule Radiation 
i "= ~~ 
> Intensity of one emission line: - 
N'A. ,. AE. ,. l6m°cv" = S. : 
e= > > — N' R 7... . ! J"A \ 
An 3 ee) a 2J'+1 
| ee, 
Particle density Electronic \. Hon! Léndon Factor 
in the levele’V’ J’ transition moment ~*~ _ 


7 
os 


=< 


¢ There are (2J+1) allowed rotational states 
(rotational multiplicity) 


¢ The selection rules allow for AJ=-1,0,1 
¢ The H6nl London Factor controls 


how the (2J+1) states are distributed 
among the branches 
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6 Emission of Atoms and Molecules tk 


Line Broadening 
Emission and absorption lines can be broadened through different processes: 


¢ natural broadening (Heisenberg uncertainty principle) 
¢ Doppler broadening (thermal motion 
* Collision broadening 


- Van-der-Waals: collisions with neutrals 
Resonance: collisions with the ‘like’ species particles (perturber’s state connected by 


an allowed transition to the upper or lower state of the transition under consideration) 
- Stark for interactions with ions or electrons, specific for each transition/line 


¢ (Zeeman effect - interaction with magnetic fields) 


fie rer enamel IN 
¢ The integration over many photons with 
i | tt pi | ty 


different shifts produces a broadening. 
> distributing the transition energy over a 
certain wavelength/wavenumber range 


¢ Temperature measurement 
(low pressures, high temperatures) 


¢ Electron density measurement 
(lines with large Stark coefficients, e.g. H, or Hs) 


2 0 2 4 6 8 10 


* An=no-N 86 
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Line Broadening 
Emission and absorption lines can be broadened through different processes: 


¢ natural broadening (Heisenberg uncertainty principle) Lorentz 
¢ Doppler broadening (thermal motion Gauss 


* Collision broadening Lorentz 
- Van-der-Waals: collisions with neutrals 
- Resonance: collisions with the ‘like’ species particles (perturber’s state connected by 
an allowed transition to the upper or lower state of the transition being studied 


¢ (Zeeman effect - interaction with magnetic fields) 
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Intro to Absorption Spectroscopy in 
Nonequilibrium Hypersonic Flows 


Canberra 


Craig Johansen, The University of Calgary, Canada 


Sean O’Byrne, University of New South Wales 
Canberra, Canberra, Australia 


sensor output (volts) 


a2 1a 
#o3 7604 7605 7606 7607 7608 7609 761 7611 
Wavelength (nm) 
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Tunable Diode Laser Absorption Spectroscopy ey 


AUSTRALIA 


Canberra 


Absorption is described bv the Beer-Lambert Law: 
—dl,, = l,k,(v)dx 
where k.(v) is the spectral absorption coefficient. 


This leads to exponential decay of laser intensity: 
Lyx — yoo * 

Measure J, know x, directly measure k,(v) 

k,(v) contains information about temperature, 


concentration, velocity of gas 


| TDLAS: Strengths and Weaknesses is 


Strengths 


Can probe non-luminous flows 
and ground-state populations 
Simple, rugged optical 
arrangements with no moving 
parts 


Comparatively inexpensive 


Fast wavelength scanning 
capability (wavelength agile) 
Can provide quantitative multi- 
parameter measurement with a 
single instrument 
- No systematic errors due to 
quenching 


Canberra 


Weaknesses 


Path-integrated technique 


— Distributions can be 
determined tomographically, 
but time consuming and 
analysis can be complex 


Spatial and temporal flow 
nonuniformity can produce bias 
errors in measurements unless 
the freestream is well 
characterized 

Some important diatomics (eg N,) 
do not have a dipole moment and 
do not have a direct absorption 
spectrum 
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| Broadening Mechanisms ux 
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¢ Spectral shape of k,( v) is from “broadening” 
¢ Pressure (Homogeneous) Broadening 
_ Avy I 
HO (VV) + (AV, 2) 


— Caused by collisions with other atoms/molecules, 
a high-pressure effect causing a Lorentzian shape 
— Generates a broadening | 
. P is the pressure 
Av, = PLL y) 7 \s the concentration 
7 was y,is collision width per unit pressure 
and a shift in the transition 6;\s collision shift per unit pressure 
Av, =P>(7,6;) 
lf temperature can be independently determined, 
pressure can be measured from these quantities 
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| Broadening Mechanisms we 
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¢ Doppler (Inhomogeneous) Broadening 
— Due to thermal motion of absorbing particles 
towards or away from laser 


Wn? —4(In 2)(v—Vo ) 
2 = Av? D | 
py m?Av, 


— Avp is the full-width at half-maximum of the 
measured peak, and is given by 


(ee 2v, yd ie —In 2 Yo is the transition center 
aan (MW )c? k is the Boltzmann constant 


| Voigt Profile ww OS 
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¢ The Voigt profile is a convolution of the two 
previous broadening effects 
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é Laser Sources vO 
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¢ In 1970s, dye lasers were often used in visible and NIR 
absorption studies, and Pb salt for IR. Now mostly 
replaced with: 


¢ Diode laser sources 

— Distributed feedback (DFB) lasers 
¢ Tune over “0.1 nm 
¢ Robust and well tested in telecoms 
¢ Very spectrally narrow (1-30 MHz) 

— Vertical cavity surface-emitting lasers (VCSELs) 
¢ Tune over 1 nm 
¢ Tune very rapidly 
¢ Spectrally narrow (<30 MHz) 

— Quantum cascade lasers (QC) 
¢ High-power pulsed source operating in the IR 


| Detection Schemes 


Canberra 


¢ Direct absorption 
— Simple, direct measurements, but poor use of dynamic 
range 
¢ Difference amplification 
— Removes intensity modulation, but only when perfectly 
balanced 
¢ Log-ratio detection 
— Simple circuit, and cancels common-mode electronic and 
laser noise very efficiently 
— Can be autobalanced, to remove offset caused by 


mismatch in reference and signal intensities 
¢ M.G. Allen et al. Applied Optics 34.18 (1995): 3240-3249 
— Bandwidth decreases as photocurrent decreases 


AUSTRALIA 


| Detection Schemes www 


Canberra 


¢ Wavelength modulation 

— Removes DC offset from signal 

— Signal proportional to the derivative of the input intensity 

— Phase-sensitive detection at a multiple of the modulation 
frequency can reduce 1/f noise in the signal 

— Traditionally required careful calibration to samples under 
controlled conditions, but recent innovation using the ratio 
of second-harmonic (2f) and first-harmonic (1f) signals to 


account for the laser intensity 
° GB Rieker et al. Appl. Opt. 48.29 (2009): 5546-5560. 


¢ Cavity-based methods 
— Useful for measurements of trace quantities, but practical 
challenges in hostile environments. 


Detection Schemes UNS 
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Direct absorption 


<e di/dA 
: 1f WMS 
j\ 


signal PD 


lock-in amplifier 


d?I/dA? 
~ 2f WMS 


current ramp 


reference PD lref - Isig 


—_—_—_—— 4 ky difference 
——— ; ; 
Sy > difference amplifier 
(+) laser 
t 
modulation - 
log (Iref / Isig) 
absorbance 
: log-ratio amplifier 
A 


98 


Detection Schemes JZ, 
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¢ £.g. Goldenstein et al. 
— Measured spectra of water vapor in UVa scramjet test facility 
— Small amplitude of second feature limits signal-to-noise ratio for direct 
absorption 
— 2f/1f ratio measurement provided consistently better signal-to-noise ratios, 
but both methods consistent with computational predictions 


O12 0.14; 
Experimental absorbance 
0.1} | ——Bestfit Voigt profile i 
0.08} 3 
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& 0,06} A 
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C.S. Goldenstein, |.A. Schultz, J.B. Jeffries and R.K. Hanson “Tunable Diode Laser Absorption Sensor for Measurements 
of Temperature and Water Concentration in Supersonic Flows,” 49'" AIAA Aerospace Sciences Meeting including the New 
Horizons Forum and Aerospace Exposition, Orlando FL, Jan 4—7, 2011, AIAA Paper 2011-1094. 


99 


Operation 
Mode 


11% Steam 


(Combustor 
Entrance) 


9% Steam 
(Exit Plane) 
12% Steam 
(Exit Plane) 
H,-Air 
Combustion 
gy = 0.33 


(Exit Plane) 


Expected Value 


700-1000 K 
11.4 + 0.2% HO 


700-1000 K 
9 + 0.2% H20 
700-1000 K 
12 + 0.2% H20 
1800-2200 K 


13% H20 


Detection Schemes une 


DA 


776 + 10K 
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860 + 30K 
9.1 + 0.2% H20 
875 + 50K 
12.1 + 0.5% H20 
1802 + 94K 
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WMS 


742 +9K 
10.8 + 0.1% H,O 


831 +9K 
9.1 + 0.1% H20 
850 + 6K 
11.5 + 0.1% H20 
1765 + 41K 


13.3 + 0.3% H20 


C.S. Goldenstein, |.A. Schultz, J.B. Jeffries and R.K. Hanson “Tunable Diode Laser Absorption Sensor for Measurements 
of Temperature and Water Concentration in Supersonic Flows,” 49'" AIAA Aerospace Sciences Meeting including the New 
Horizons Forum and Aerospace Exposition, Orlando FL, Jan 4—7, 2011, AIAA Paper 2011-1094. 
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| Velocity Measurements 


¢ Simplest quantity to measure in a high-speed 
flow 
— Measured through Doppler shift between beams 
with velocity components in two different 
directions relative to the flow 


We CAV no pp 


V Soupee 
— Shift in the peak position relatively insensitive to 
amplitude noise. 
— Shift typically calibrated to etalon measurements 
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| Velocity Measurements — wsw 
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¢ F.g. Lyle et al., mass flux sensor based on measurements of 
density and velocity using oxygen A-band near 760-nm. 
¢ Measured speed with time-averaged precision of 0.25 m/s. 


1.00 
> 0.99 
j. 
oO 
5 0.98 
€ 
Boamn 4 z 0.97 Beam 1 
Air flow \a a ——-— Beam2 


Velocity = U A 
20 


/ \ 


DAQ computer 


13083.6 13084.0 13084.4 13084.8 
Frequency, v [em] 


K. H. Lyle, J. B. Jeffries and R. K. Hanson “Diode-laser Sensor for Air-mass flux 1: Design and Wind Tunnel 
Validation,” AIAA Journal, 45 (9), 2204-2212, 2007. 
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| Translational Temperature 
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¢ Can be determined through the Doppler width, but 
AS Trans © a“e , small errors in Doppler width have a 
significant effect on T 
AT  2A(Av,) 
T AV, 

— F.g. An uncertainty in width measurement of 1.5% for a 
Doppler-broadened transition in the oxygen A-band will 
cause an uncertainty of 3% in temperature 

— Uncertainty further increased in the presence of pressure 
broadening 

— Measuring T through linewidth has the advantage of only 
requiring 1 transition for a measurement 
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| Translational Temperature 
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¢ F.g. Measurements of 
temperature in an . Seg 
arcjet facility using an i. / —_a 
AR I line at 811.531 nm 4 = | | 

¢ 12% uncertainty in ° 
width translated to 30% 
uncertainty in 7,,,,, over 
2000—1000 K 


rans 


Time (ms) 


F-Y Zhang, K Komurasaki, T lida, and T Fujiwara, “Diagnostics of an Argon Arcjet Plume with a Diode Laser,” 
Applied Optics, 38 (9), pp. 1814-1822, 1999. 
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é Translational Temperature 


1 COC BW 


¢ F.g. Measurements of 
NO translational 
temperature at 5.44 um 
in the LENS | facility, 
using direct absorption 


significantly different to 
the predicted 
temperature of 563K 


Fit temperature of 290 K 
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& TOO 


—Theory 290K 
© Data 


1835.58 


1835.56 


1835.57 
Wavenumber (cm-1) 


AOC 
fe) 
° 


— Theory 563K 
© Data 


1835.58 


R. A. Parker, T. Wakeman, M. MacLean, and M. Holden, “Measuring Nitric Oxide Freestream Concentration using 
Quantum Cascade Lasers at CUBRC,” 44" AIAA Aerospace Sciences Meeting and Exhibit, AIAA Paper 2006- 


926, 2006. 
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¢ Typically measured using the ratio of integrated 
absorbance a over two lines with different 
linestrength variations with T 
a,=Py7LS, 
— for transition j of an absorbing species / 


— Lis the path length and Sis the linestrength of the 
transition 


¢ All quantities other than S, are common to the two 
lines and divide out 
hcAE,,,( 1 1 
ay = S; _ Si) | k ie 4 
a, S, S,(Z) 


R. A. Parker, T. Wakeman, M. MacLean, and M. Holden, “Measuring Nitric Oxide Freestream Concentration using 
Quantum Cascade Lasers at CUBRC,” 44'" AIAA Aerospace Sciences Meeting and Exhibit, AIAA Paper 2006- 
926, 2006. 
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Rotational Temperature 


Canberra 


E.g. Wehe et al. measured H,O 
rotational and translational 50x10" 
temperature in the 10 MJ/kg 40 
Calspan 96-inch hypersonic shock 
tunnel freestream 

v1 + v3 band near 1396 nm (A,) 


30 


Absorbance 
Nm 
oO 


and 1400 nm (A,) "0 
8 kHz scan rate 0 eee 
Tot uncertainty of +2.3 % 0.00 0.05 0.10 0.15 0.20 0.25 


. . . “41 
compared with 7,,.,. uncertainty Relative frequency (cm ) 


of 2.7 % and +6.4 % for the T,,..,,. 
measurements on the two lines 
S. Wehe, D. S. Baer, and R. K. Hanson “Tunable Diode-laser Absorption Measurements of Temperature, Velocity, 


and H.O in Hypervelocity Flows,” 33 Joint Propulsion Conference and Exhibit, Seattle, WA, AIAA Paper 97- 


3267, 1997. 
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Vibrational Temperature 


2100 pe 


E.g. Brandt and Roth measured T,,, 
and T,,,,. in CO in a reflected shock = 
tube, using a lead salt laser operating on 
at 2100-2200 cmt, at 20 kHz scan 


rate. 


vibrational temperature Ty 
i 
KT al 


— calculation 


T4as 
a é 7 0 400 800 1200 1600 1S 2000 
— aan time — time particle time 
O. Brandt, and P. Roth, “Temperature Measurements behind Shock Waves using a Rapid Scanning IR-diode 
Laser,” Physics of Fluids, Vol. 30, No. 5, pp. 1294-1298., 1987. 
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TDLAS Tomography 
Major limitation of TDLAS: path-averaged 
— Can do tomography 


36.3 mm 
et Lo 
od 36.3 mm 


_ t 


Weta! 


ebdbdds 


coaxial to 
23 24 25 26 27 28 29 30 DAD 


— 15x15 grid provides data at 225 points (36 mm grid) 
L. Ma, X. Li, S. T. Sanders, A. W. Caswell, S. Roy, D. H. Plemmons, and J. R. Gord, “50-kHz-rate 2D imaging of 
temperature and H,O concentration at the exhaust plane of a J85 engine using hyperspectral tomography. 
Optics Express, Vol. 21, Issue 1, pp. 1152-1162, 2013 
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| TDLAS Tomography 


¢ General Electric J85 gas turbine engine 


— University of Tennessee Space Institute (UTSI) 
co 


Canberra 


— 60m ethernet 


TDM 3-FDML 


—60m BNC ______ 


4x 32 A rl. 
=60m 1 — fF Taw wl exall 


Control Room ot 


e Measure water and temperature at 50 kHz rate 


L. Ma, X. Li, S. T. Sanders, A. W. Caswell, S. Roy, D. H. Plemmons, and J. R. Gord, “50-kHz-rate 2D imaging 
of temperature and H.O concentration at the exhaust plane of a J85 engine using hyperspectral tomography,” 
Optics Express, Vol. 21, Issue 1, pp. 1152-1162, 2013 110 
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| TDLAS Tomography Results he: 
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¢ Measurement is over ~0.5 meter square 
— Time resolution is 20 usec (50 kHz) 
— Spatial resolution is ~36 mm 


L. Ma, X. Li, S. T. Sanders, A. W. Caswell, S. Roy, D. H. Plemmons, and J. R. Gord, “50-kHz-rate 2D imaging of temperature and H,0 7; 4; 
concentration at the exhaust plane of a J85 engine using hyperspectral tomography,” Optics Express, Vol. 21, Issue 1, pp. 1152-1162, 2013 
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L. Ma, X. Li, S. T. Sanders, A. W. Caswell, S. Roy, D. H. Plemmons, and J. R. Gord, “50-kKHz-rate 2D imaging of 
temperature and H,O concentration at the exhaust plane of a J85 engine using hyperspectral tomography,” Optics 
Express, Vol. 21, Issue 1, pp. 1152-1162, 2013 112 
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| Using Spectral Information ee 
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¢ Many hypersonic environments (flight vehicles, test 
facilities) have limited optical access for tomography 
— Can use spectral information to infer spatial 
variations of temperature. 


sy" 296, 712 K 


T=496K | 
, | 


absorbance per cm 


13140 13145 13150 13155 
“4 
frequency, [cm ] 


S.T. Sanders, J. Wang, J.B. Jeffries, and R. K. Hanson. "Diode-laser absorption sensor for line-of-sight gas 
temperature distributions." Applied Optics 40,(24) 4404-4415., 2001. 
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Extracting Line-of-Sight Temperatures ,, sy 
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e Example: Temperature in a tube furnace using a single VCSEL line of sight 


e The ‘distribution technique’ assumes a priori knowledge of the distribution 
shape. 


Temperature 


Absorbance 


Distance | J * 


Dal = 2 ee ee a a a rik] Wavenumber 


‘experimental’ spectrum ——— 
theoretical spectrum — — - 


L. C. Webster, S. O’Byrne, A. F. P. Houwing, “Determination of Temperature Distributions in Air using a 
Scanning Vertical Cavity Surface-Emitting Laser,” Proc. 4° Australian Conference on Laser Diagnostics in 
Fluid Mechanics and Combustion, Adelaide, Australia, 2005. 114 
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| TDLAS Hyperspectral Tomography usw GY 
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¢ Using hyperspectral sources, can gain spectral images 
over many lines, reducing number of projections 


aon Spectra Measured at Beam Location 6 
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L. Ma, W. Cai, A. W. Caswell, T. Kraetschmer, S. T. Sanders, S. Roy, and J. R. Gord, “Tomographic imaging of 
temperature and chemical species based on hyperspectral absorption spectroscopy,” Optics Express 17(10), 8602— 
8613, 2009. 115 


Flight-Testing TDLAS Systems ww @Y 
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¢ Small size and power requirements have allowed 
TDLAS systems to be used in hypersonic flight tests 


T. Barhorst and S. Williams, "Development of an In-Flight Non-Intrusive Mass Capture System," 45% 
AIAA/ASME/SAE/ASEE Joint Propulsion Conference & Exhibit. AIAA, Denver, Colorado, USA, 2009. 

J. Kurtz, M. Aizengendler, Y. Krishna, and S. O’Byrne, S (2015) “Flight Test of a Rugged Scramijet-Inlet 
Temperature and Velocity Sensor,” AIAA SciTech, 5-9 January, Kissimee, Florida, USA. 116 


| Flight-Testing TDLAS Systems usw 


- 


Canberra 


Relative Signal 


250 


200 


150 


100 


Wavelength Modulation Absorption 


+ signal 
| O reference 


sensor output (volts) 


6. Py si io id id ia id t r 
10 20 30 40 50 60 70 %o3 760.4 760.5 760.6 760.7 760.8 760.9 761 761.1 
Point Number Wavelength (nm) 


M. S. Brown,and T.F. Barhorst. "Post-flight analysis of the diode-laser-based mass capture experiment onboard 
HIFiIRE flight 1." In Proceedings of 17th AIAA International Space Planes and Hypersonic Systems and 
Technologies Conference paper AIAA-2011-2359. 2011. 
J. Kurtz, M. Aizengendler, Y. Krishna, and S. O’Byrne, S (2015) “Flight Test of a Rugged Scramjet-Inlet 
Temperature and Velocity Sensor,” AIAA SciTech, 5-9 January, Kissimee, Florida, USA. 
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TDLAS is a reliable and well established nonintrusive 
measurement technique, capable of precise path- 
averaged measurements of velocity, density and 
trans/rot/vib temperatures 

Availability of telecoms lasers and detectors makes it 
comparatively economical 

Path-averaging disadvantages can be mitigated using 
tomography or additional spectral information 
Particularly valuable for flight-test applications and 
time-resolved online facility monitoring 

Rapid development in mid-IR sources and detectors 
will significantly expand applications 
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Intro to Planar Laser-Induced Fluorescence for ee 
Hypersonic Nonequilibrium Flows ware 


Brett Bathel, NASA Langley Research Center, Virginia, USA 
Paul Danehy, NASA Langley Research Center, Virginia, USA 
Craig Johansen, The University of Calgary, Canada 
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Figures from Inman, J., Bathel, B., Johansen, C., Danehy, P., Jones, S., Gragg, J., and Splinter, S., 
(2013) “Nitric oxide planar laser-induced fluorescence measurements in the Hypersonic Materials 
Environment Test System,” A/AA Journal, Vol. 51, No. 10, pp. 2365-2379 
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| Planar Laser- Induced Fluorescence (PLIF) “” ® 


“Tunable Laser 


Laser sheet 
excites molecules 


Excited molecules 


fluoresce 
Excited state 


CCD 
camera 


Ground state 
detects 


LIF ~ nyo 
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é PLIF can probe nonequilibrium & & 
S.eeees—t iF 


Commonly used in equilibrium environments 
Useful probe for thermal non-equilibrium: 
— Translational temperature 
— Rotational temperature 
— Vibrational temperature (often 7,,, 4 7,04 = Terran) 
— Species specific: 
¢ Different species can have different temperatures 
(especially vibrational temperatures) 
¢ Can measure, quantify species in chemical 
non-equilibrium (eg. NO in shock tunnels) 


| LIF Theory (Two-Level Model) © Gy 
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D Energy of absorbing species at excited state (2) 


Energy of absorbing species at ground state (1) 


Assumptions: 


A. There are only two important states to consider (ground state 
and excited state). Other levels (rotational and vibrational) are 
ignored. 


B. Prior to excitation, population of absorbers in excited state is 
negligible (N,=0) 
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| Two-Level Model Theory & & 
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W,.= Stimulated Absorption 


An incident photon from the light source (e.g. laser) is absorbed and 
sends some of the ground state population to the excited state 


dN, 
—— = —N,\W,, : 
dt Population rate equations 
incomplete at this point 
dN, ( p point) 
= NW, 


dt 123 


| Two-Level Model Theory & & 


Q,.= Collisional Excitation 


The excited state population is further increased by collisional 
excitation. 


dN, 
We = —N,(W,, + Qi2) 


dN, 


dt. — NW, + Q12) 
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| Two-Level Model Theory ; GY 


W,,= Stimulated Emission 


Some of the excited state population returns to the ground state 
through stimulated emission 


dN, 
dt 


dN, 
“dt. = NW, — Qi2) — NW, i 


= —N,(W,,4+ Q12) + NW, 


| Two-Level Model Theory & & 


Q,,= Collisional Quenching 


Some of the excited state population returns to the ground state 
through collisions (self or with other atoms/molecules) 


dN 


777 = —Ni{(W2 + Qi2) + N2(Wa1 + Qa) 


dN, 


Fe = NW + Q12) — NoWar + Qa1) ts 


| Two-Level Model Theory & & 


Q,.,e= Predissociation 
W on= photo- 
ionization 


The excited state population gets depleted through dissociation and 
ionization (without returning to the ground state) 


dN, 
‘dt —N, (Wi y2 + Qi2) + NW; + Qa1) 


dt — N,(W, + Q12) i NW, + Qo + nr + W ion) 
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| Two-Level Model Theory & & 


Qhre 


ion 


A, 


A>,= spontaneous 
emission 


Finally, some of the excited population return to the ground state via 
spontaneous emission and emit a photon (fluorescence) 


aN 
a = —N,(W,2 + Qi2) + N2(W 2, + Qo, + Any) 


aN 
= Ny (Wy + Qi2) — No(Wo1 + Qa1 + Qpre + Wion + Ary) 
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é Two-Level Model Theory , & & 
— Ke 


Fortunately, for many applications, Q;, Ws, Qi.,, and Q,,. are negligible. 


dN, 
aa, oo = NW, + Qo, + Az) — N,(W,, +42) 


dN 


o - N,W,, +Q%2) a NW, + Qo1 t Are EWion = Ai) 129 


| Two-Level Model Theory & & 
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sircliestion cies fon: @Ni _ _ aN 
implification of rate equation: rs — 7s = N(Q>, + A>1) —N,W,, 
Mole fraction of absorbing species Total number 
, density 
Conservation of mass for N,+N,=N.=x.fpN 
absorbing species: ' 2 ° ° atte 
Number density of the absorbing Boltzmann fraction 


species in the ground state prior to excitation 130 


é Steady State Two-Level Model @y 


CALGARY 
Steady state assumption aN, _ aN, = 
(e.g. continuous wave laser): dt dt 
Xsf pNW12 


Solve for N.: Nz = W12+0>,+421 


Weak excitation assumption: W,, « (Q,, + A>) 


Caet 
Integrate spontaneous emission | N,A,,dt = NoAo tae, 
rate: 0 


volume 


ayeiied solid angle 


Include the optics of detection detection 
to get LIF signal: time 


Sue = Xef gNW — )t Vy 
Pa LIF sJ BY“ T"" 12 det 4 


total number of 7 ~4i 21 collection 
photons collected fluorescence yield (®) efficiency 


é Concentration of Absorbing Species & & 
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0) 
Sip = XsNrf pW 2P Cae FN 


Figure from Inman et al (2013) 


¢ To first order, S,,- is proportional to the 
concentration of the absorbing species 


¢ Qualitative PLIF visualization often used 
to show the spatial distribution of 
absorbing species 


¢ (see the text for measurement strategies Sir X& XnoNr 
for quantitative concentration) 
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é Boltzmann Fraction & & 
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0, 
Sip = XsNrf pW y2PtacV Aq 


Rotational Boltzmann fraction 
Computed using LIFBASE 


f ge 5B (J. Luque, SRI International) 
CW Ej/ kT 
2 di hd i/ B 

| | 100 K 


— a function of temperature only I, 
— Peak at 300 K is ~3% for N = 7 : 
(J = 7.5) 
¢ Vibrational and electronic 
Boltzmann fractions = = —sARTEEEE ttn... 


— Most molecules in ground state for 1000 K 
T < 2000 for atoms, diatomics alt 


Population 
— 
ee) 
oO 
=) 
A 


Rotational quantum number, N (or J) ;33 


©) 


Stimulated Absorption & & 
: CALGARY 


() 
Sur = XsNrf BW 1 2PlgeV An 


W,. is the probability of an absorption transition per sec 


absorption 
lineshape 
laser irradiance function laser spectral 
profile 


W,,= Bp! | Y,L.dv 
A 


Einstein B coefficient Y 


Pa 


spectral overlap integral 
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| Fluorescence Yield 


UNIVERSITY OF 
CALGARY 


C), 
Sip = XNrf pByIGPtz.V ae 


The fluorescence yield is affected by collisional quenching 


A., —— Probability of an atom/molecule 
om = 21 undergoing spontaneous 
Q>, a2 Ass, emission per second 
Probability of a quenching ~ (Einstein A coefficient) 
collision per second 


Collisional quenching is hard to predict since it depends on 


many factors wt 
Nr a P/kpT 


Qo = Nr dLuiX Fs isi y Is mole fraction 


o iS quenching cross section 
v is mean relative velocity 


135 


| Steady State S,,, dependencies & & 


Str % XsfeMT)By2IGYs5,P,T, VPNs, Xi, PT )taet 


— LIF signal is sensitive to 7} PB U, 7., 7, 
¢ Can potentially measure these parameters 
¢ Complicates measurement of a single parameter 


— Making quantitative measurements with LIF 


involves finding a way to make measurement 
sensitive only to the parameter being measured 
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é Time Dependent S,- Solution ® & 
9 ilime ve CALGARY 


S7 7A) = Sr (Cblaserde ~ (ttlaser)/TLIF t > bias 


x 
photons oe steady ; 
collected Trp = (Qo, + Agi)” 
per time 
exponential 
_-decay after 


laser pulse 


Stir) 


t 


— Laser pulse length is typically ~10 ns 
— LIF lifetimes can be a few ns to 100’s of ns or longer 
¢ If you can measure lifetime, you can measure Q! 
¢ Can measure velocity from long lived decay (MTV) 
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| More Complicated LIF Models & & 
CALGARY 


¢ More complete LIF theories and phenomena: 
— Saturation effects (see manuscript) 


— Line shapes 

¢ overlap integral, G, to be discussed below 
— Multi-level LIF models 2 

¢ 2 photon absorption i ial 

e 3,4, 5 level models i 


— Predissociation 

— Rotational energy transfer 

— Vibrational energy transfer 
¢ For further reference: 


A. C. Eckbreth, Laser Diagnostics for Combustion Temperature 
and Species 2"¢ Ed., Combustion Science & Technology Book 
Series, Volume 3, Taylor & Francis, New York, 1996 


1 
2-photon LIF, 3 level model 
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Arealmolecule:NO  @& 


[10°cm |] 
¢ NO is a diatomic a Fas. 
molecule ' ee 
— Electronic ~ 
— Vibration 40 
— Rotation ; 
¢ Not shown electronic 45 
¢ Absorption and “2° . ; 
fluorescence take 20 ~n 
place between — nile 
electronic states " 
— Show rotational and ; Ws 
vibrational structure . 15 20 Al 


Internuclear Distance 


Figure A2.1. Selected potential energy curves that show the electronic structure of NO. In this 
experiment. we probed transitions in the 4A—Y band of NO [adapted from (Scheingraber and Vidal.1985 


LIFBASE Simulation of NO 


UNIVERSITY OF 


CALGARY 


“Bi urease 2.1.1 Sen) 
—— Oo JT a oe WO &-O-S-@8-?-B-y ~F- AVHRR 


Database Simulation | LIF/Emission settings | Other options | Display | Experimental data | 


; Range in A (Vacuum) Instrumental 4 (¥acuum] Total resolution (instrumental + broadening) = 0.0427 A (Vacuum) 
Kind | Absorption _¥ Start — |2265.5 Resolution | 0.0221 \¥ Doppler broadening = 0.00512 |¥ Predissociation broadening 


Units |A (vacuum) + End 2270.5 Lineshape | Gaussian 7] IV Collisional broadening |0.6 = *{P/101.32)"(300/T)"|1 = 0.03063 
Noise (%] |0 Pressure |101 kPa | Lorentzian(%) [0 Temperature | 300 K |v Thermalized system Population distributions 


Simulation NO (A-X) Absorption 


LIFBASE 
freeware by 
J. Luque, SRI 
International 


6. ihe he: He neta lifbase@yahoo.com 
4 6 6 TF [8 igo 


4 13 2 Wi fio 3 is vi 


b12) 


20 1 
O12(2) 13 


7 |6i5 
Ai 212) 


2266.0 2266.5 2267.0 2267.5 2268.0 2268.5 2269.0 2269.5 2270.0 
x= 2268.700, y= 1.513 Wavelength (A) (Vacuum) 


LIF Excitation Scan between 226.6-227 nm, 1 atm, 300 K 
— Lines all have band names (O, P,Q, R, S) and N’s or J’s (rot quant #) 140 


| LIF Temperature Measurements © @y 
cr 


¢ Rotational Temperature 


— Ratio of 2 or more rotational lines 
— Describe general case 


— Example in a turbulent base flow 


¢ Vibrational Temperature 
— Ratio of 2 or more vibrational lines 


— Example of freestream characterization 
¢ Boltzmann plot, where more than 2 lines used 


¢ Translational Temperature 


@ Rotational Lines of NO, 300 K 


Duras aii 


—<—r GF 


UNIVERSITY OF 


CALGARY 


7a isoex OO |e-O-S-B8-?- By \~F AVE 


Database Simulation | LIF/Emission settings | Other options| Display | Experimental data | 


Range in A (Vacuum) — Instrumental A (Vacuum) Total resolution (instrumental + broadening) = 0.0315 A (Vacuum) 
Kind Absorption +] | stat [2255 Resolution {0.001 I Doppler broadening = 0.00512 I Predissociation broadening 
Units |A(Vacuum) v] End = [2277 Lineshape [Gaussian | || [¥ Colisional broadening /a.g P/101.32p(300/T}"[1 = 0.03063 


Noise (%) [0 Pressure [101 [kPa] | Lorentzian (%) [0 


Temperature[300  K [¥ Thermalized system Population distributions 


Simulation NO (A-X) Absorption 
I 


Vy 


I 
= 2259.802, y= 62.168 


Wavelength (A) (Vacuum) 


e LIF Excitation Scan between 225-227 nm, 1 atm, 300 K 
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@ Rotational Lines of NO, 600 K & 


Eh urease 2.1.1 She aes 
B: Sa |noow IO | ae- Oy S-@B-?~-B-|~rF AVA 
Database Simulation | LIF/Emission settings | Other options| Display | Experimental data | 
Range in A (Vacuum) — ~ Instrumental A (Vacuum) Total resolution {instrumental + broadening) = 0.0321 A (Vacuum) 
¥ Doppler broadening = 0.00725 I Piedissociation broadening 


Resolution | 0.0221 


Kind Absorption ~] | stat [2255 
Units [A (vacuum) +] End — [2271 Lineshape [Gaussian +] | M Colisional broadening [0.6 {P/101.32"(300/T/"[7 = 0.01531 
0 Temperature[60) kK [¥ Thetnalized system Population distributions 


Lorentzian (%] 


Noise (2) [0 Pressure [101 [kPa >] 


" Simulation NO (A-X) Absorption 
| 


2270.0 


et LP PAY UU UY hae ©, 


2260.0 
Wavelenath (4) (Vacuum) 


| | 
2266.0 


x= 2258.116, y= 105.000 


e LIF Excitation Scan between 225-227 nm, 1 atm, 600 K 
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é LIF: Rotational Temperature, 7 a) 
; CALGARY 


YO 
¢ Recall: 


Sup © XsfaT)By2IGQ5,P,T, VOY Xi PT )taet 
¢ Obtain LIF or PLIF on two different rotational 
lines of a molecule: i and 7. Take ratio, R: 


SLIFi __ R=C By jE; (2Ji+Dexp[-Fy i/kpTrot] 


SLIFj By2 jE; (2Jj+1)exp[-Fy j/KpTrot| 
¢ Where: 
QO, x. 1. G, V; 
— Cis aconstant N,and ® 
_ ZF: . (including Q,,) 
E is the laser energy (in mJoules) Ai Res 


— J is the rotational quantum number 
— F, is the energy of the state 
— kz is Boltzmann’s constant 144 


YO 


| LIF: Rotational Temperature, 7 oO ay 
CALGARY 


¢ Solve for 7. 


rot’ 


—AErot 


palin) 
kpln| CRD 
: | By2 jEj(2Ji+1) 


(ee — 


¢ Where: 
AE rot = (Fy — Fj) 
¢ Temperature sensitivity: 


OTrot _ KpTrot OR 
Trot AErot R 


¢ Suggests to maximize AE,,.; (as long as signals still 
strong) to minimize error in temperature, oT 


rot 


— Should choose widely separated states: e.g. /=3.5&35.5 |. 


é Typical Strategies for 7, & & 


¢ Typically two approaches: 


— Steady laminar flow 
¢ Obtain images on one J, then scan laser to another J 
¢ Ratio images, determine temperature 
— Turbulent flow 
¢ Two lasers, two cameras (Cumbersome! Expensive!) 
¢ Acquire two images nearly simultaneously, different J’s 
¢ Need to do this because ratio is a nonlinear function of 
temperature: Suri _ p — p B12ibi (2Jit Dexp|-Fyi/kpTrot] 
Supj «Baa, jE (2Jj+1)exp[-Fy,j/kaTrot] 
¢ Typical single-shot temperature measurement precision 
is ~20% at flame temperatures (Seitzman, 1994) 


¢ If you use 1 laser/1 camera, >bias errors in turb. Flows 
(E. R. Lachney, and N. T. Clemens, “PLIF imaging of mean temperature and 
"pressure in a supersonic bluff wake,” Experiments in Fluids, 24, pp. 354-363, 1998.) 
— Average many measurements to reduce uncertainty 
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@ sinale laser/camera for T.,,in turbulent flow © Gy 
CALGARY 


¢ Lachney and Clemens (1998) used 1 laser / 
1 camera to map mean 7, in a base flow. 

¢ NO seeded N, at Mach 3 

¢ Chose: J = 8.5 and J = 10.5 


M=3 


2.0 pe 
18 Recompression 
) Boundary Expansion Sneek 
1.6 layer fan 
1.4 “\ 
12 shee 
WN ee 
© 1.0 2 ~ re = 
0.8 i = — 
0.6 
0.4 Splitter plate Recirculation 
0.2 zone Turbul 
: urbulent 
50 100 150 200 250 300 
Temperature (K) Turbulent supersonic wake flow 


Ratio of lines, used to measure 7,,, ‘a 


La 


Ge 


, > ~ ate fg ~ ft 
©, Single laser/camera for T,,,in turbulent flow“ (as 


<a 


Lachney and Clemens (1998): 


Flow 


Turbulent supersonic wake flow temperature measurement 
Agreed with temperatures inferred from probe measurements to within 11% 


E. R. Lachney, and N. T. Clemens, “PLIF imaging of mean temperature and 
pressure in a supersonic bluff wake,” Experiments in Fluids, 24, pp. 354-363, 1998 La 


Vibrational Bands of NO 


——— Oo OT aR [soon 10 |\e-O-~ S-O-?-B- Wr AVHRR 
Database Simulation | LIF/Emission settings | Other options | Display | Experimental data | 

: Range in (Vacuum) - Instrumental A (Vacuum) -- Total resolution {instrumental + broadening) = 0.0296 A (Vacuum) 

Kind Stat [2200 Resolution [0.0221 ¥ Doppler broadening = 0.00971 PM Predissociation broadening 
sian >] ||  Collsional broadening [0.6 —(P/101.32y(300/T}"[T = 0.00931 


Noise (%) [0 Pressure [101 kPa >] | Lorentzian (%) [0 Temperature|1000 kK [¥ Thermalized system Population distibutions 


Absorption _¥ | 


Units [A (vacuum) +] | End — [2500 Lineshape [Gau: 


| | I 
2350.0 
x= 2276.405, y= 103.053 Wavelenath (4) (Vacuum) 


¢ Three vibrational bands: T = 1000 K 


— Vibrational bands have quantum number, v 
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UNIVERSITY OF 


CALGARY 


Vibrational Bands of NO 


= 9. Fa Roe 10 s-O-S-@-?-B- ~-AVaRR 
Database Simulation | LIF/Emission settings | Other options | Display | Experimental data | 

Range in A (Vacuum) — - Instrumental A (Vacuum) — 
Kind Absorption 


Total resolution (instrumental + broadening) = 0.029 A (Vacuum) 
Start 2200 Resolution | 0.0221 Doppler broadening = 0.01418 


I Predissociation broadening 
 Colisional broadening [0.6 "(P/101.327(300/T)"[1 = 0.00528 


Temperature[2000 _ [¥ Thermalized system Population is! 


Units [A (Vacuum) +] | End — [2500 Lineshape [Gaussian _~| 


Noise (%) [0 Pressure [101 kPa >| Lorentzian(%) [0 


I | I 
= 2293.90, y= - 4.757 


2350.0 
Wavelength (&) [Vacuum] 


¢ Three vibrational bands: T = 2000 K 
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6) Vibrational Temperature 7), Measurement & & 
CALGARY 


¢ Asimilar approach can be made to obtain 


vibrational temperature measurements 
— Typically excite the same rotational line (J) in 2 or 
more vibrational bands (F; = constant) 
— Can do 2 line measurements as before a 
— Can use a Boltzmann Plot: ae 
¢ Graph: sates probed 
SLIF,i 


7 | versus G,,, 
B42 jEj(2Jj+Dexp|-F 7 i/kpTrot| 


¢ Slope is: lf F, = constant then 7,,, 
—1/(KpT yin) can be determined 
independent of T,., 
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Vibrational Temperature 7), Measurement & & 
CALGARY 


Vibrational temperature measured in a nozzle exit 
flow verifying vibrational nonequlibrium (T,,, 4 T,,,) 


: 1000 
4 © Experimental = 
ae Fit to data : 785 < 
. 2 ~~ — T= 850 2 900 | 
S : us 2 T yi» (LIF) = 785 +/- 30 K 
on ® g00 
ae -4 Ss 
kd o 
eg 700 
a 5 = 
LE 7 % E00 T yi, (CFD) = 670 K 
= 6 
8 > 
9 500 


260 270 280 290 300 310 320 330 


0 a oe a a Distance from nozzle throat (mm) 


Vibrational Energy, G ,» (cm’') 
T= 400 K and decreasing with distance 


P. C. Palma, P. M. Danehy, and A. F. P. Houwing, “Fluorescence Imaging of Rotational and 
Vibrational Temperature in Shock- Tunnel Nozzle Flow,” AlAA Journal, 41(9), September 2003 
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é Translational Temperature & & 
¢ Can measure T,. (and P) from width and 


shapes of spectral lines 


¢ Spectral lines have a shape 
— Molecular motion leads to: : 


: : Gap (min2)'” 
¢ Doppler Broadening (inhomogeneous) 
¢ Gaussian shape 


* FWHM: = av, =*2 pane 


— Collisions lead to: 
¢ Pressure Broadening (homogeneous) 
¢ Lorentzian shape ie 
P 295 
© FWHM: = 4ve=7T a5" 27 4,| 2 +f 
— Combination is: Voigt 
¢ Convolution (if Av,,Av, comparable) -* -: - 


Modest, M: Radiation Heat Transfer, Academic 


. Formulas In manuscript Press, Elsevier Science, USA, 2003 153 


| Translational Temperature & & 


RF 
plasma generator 


laser 


¢ If using a very narrow soa ft 
1 1 “ania | iz t , — matching 
linewidth laser... ae vt 
Av, <<Av,,AVp | : cnochromator] 
..can neglect laser width (latin ati 
— Semiconductor diode laser chai (f= ; Leer 


¢ Scan laser frequency across 


es -% @ z=1imm (r, z) = (0, 11 mm) 
a transition 2 — data 
s ehangahe 5 
— Argon atom at 810.4 nm = nye 9.2x10"%om 
— Turbulent inductively 2 
coupled plasma (ICP) torch j 
¢ Fit measured LIF spectrum 3 
— Have homogeneous (Stark) e. 
and Doppler broadening frequency [GHz 
— Determine n, and T Sirscepietis Grecaife pleaines arg semniondicior sere dined ti 


Optical Society of America B, 9(11), pp. 1968-1978, November, 1 


©) 


Species Concentrations & & 


Problem: How to deal with unknown @Q,,? 
Solutions: 
— Saturated LIF: (W,2 + W21) > (Q21 + Ao1) 
° Q,, is negligible 2 
compared to W’s 
¢ Hard to achieve 
— Predissociation 
* On7¢ >> Qo) 
— Measure 7, 1 
¢ Determine Q,, from 7, 
— Measure other species present and predict 7, 
¢ Raman/Rayleigh/LIF (Sandia and DLR) 
— Counter-propagating beams, take ratio (versluis, 1997) 


Potential problems: absorption, rad. trapping... 


| LIF Advantages/Disadvantages & & 
CALGARY 


e Advantages 
— Relatively easy to make planar measurements (PLIF) 
— Provides flow visualization simultaneously 
— Relatively easy to perform 
— Can probe many different parameters (T, V, concentration) 
¢ Velocity measurements (MTV) can have with 1-2% precision 
¢ Disadvantages 
— Precision of single-shot temperature measurements 
(~20%) is not good compared to CARS (2-3%) or LITA (1%) 
— Difficult to do quantitative concentration measurements, 
especially planar 
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6) PLIF Backup Slides 


é LIF Theory: Atoms & & 


Lyman Series 


e- 
O 
@ n=] Balmer series 
n=2 A= 656 nm 
v= 4.6x10'4 Hz 
H atom n=3 v= 15,240 cm"! 
Paschen senes 
n= 4 
n=3 n=6 
¢ Two level model Hydrogen atom energy levels 
(from Wikipedia, “hydrogen spectral series”) 
—e. g. Hyd rogen atom http://creativecommons.org/licenses/by-sa/3.0/deed.en 
¢ Photon energy =hv 
—c=hv wavenumbers =v =1/A cm? - 


Introduction to LIF Velocimetry Techniques for @@y 
Hypersonic Nonequilibrium Flows ww 


Brett Bathel, NASA Langley Research Center, Virginia, USA 
Paul Danehy, NASA Langley Research Center, Virginia, USA 
Craig Johansen, The University of Calgary, Canada 
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NASA and The Australian National University 


Texas A&M 
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é Outline: Velocimetry Techniques Ss 


¢ LIF-Based Velocimetry Techniques 

— Doppler shift 
¢ Based on laser frequency relative to absorption line 

— Molecular tagging velocimetry (MTV) 
° aka Flow tagging velocimetry 
¢ Time of flight 
e Single-laser method requires lifetime long enough for flow to move while 

light is being emitted 
— Summary 


¢ Other Techniques 
— Particle: Particle Image Velocimetry (PIV), Laser Doppler 
Velocimetry (LDV) 
— Interference/Scatting: Laser-Induced Thermal Acoustics (LITA) 
— Schlieren: Density Tagging Velocimetry (DTV) 
— Scattering: Rayleigh 
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Doppler Shift Velocimetry 


6) Doppler Shift Velocimetry Ss 


¢ The Doppler shift: 


Av =—Vo 
Cc 
Wikipedia “Doppler Shift” 


f 

* Can be used to measure velocity with LIF 
¢ Relative shift between laser and absorbing molecules 
e Angle of detection does not matter for LIF 

¢ The overlap integral, G, which describes the overlap 
of the laser’s spectrum with a line 


¢ L, is laser’s spectral shape, typically modeled as Gaussian 
or Lorentzian 


¢ Y is the transition’s absorption line shape, can be shifted 
by Doppler shift 


G = [ L,Y,dv 
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@ Doppler-based LIF Velocimetry S& 


¢ Direct laser sheet into the flow 
¢ Acquire images or signal at different 
laser frequencies 
— Choose frequencies near an 
absorption line 
— Measure absorption lineshape at 
each point 
— Reference data acquired 
simultaneously 


cylindrical 


lens 
Sy fluorescent cell 


CCD camera 
beamsplitter 


deflection 
mirrors 


Center: 0.493+/-0.03pm 226.1 39 nm| 


y [mm] 


Fluorescence intensity [-] 


oluas= = = 
15 -1 -05 0 O05 1 #15 2 25 3 35 
Laser detuning [pm] 


— Measure from two different 
directions to get two different 


velocity components R. Hruschka, S. O’Byme, and H. Kleine, “Two-component Doppler- 


7 : : shift fluorescence velocimetry applied to a generic planetary entry 
Do need to consider pressure shift probe model,” Experiments in Fluids, 48, pp. 1109-1120, 2010 
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| Doppler-based LIF Velocimetry 


Self-referenced tuning Doppler 


velocimetry 
— Does not require reference cell 
measurement 
aes ual —_¢_,__(Av/2) 
Oe cos@ ([v¢t+v;]/2) 


— Laser passed through measurement 
volume, followed by a delay line that 
returns beam along original path 

— Requires high-speed measurement 
device (e.g. photomultiplier tube) to 
resolve forward and return signals 


Fixed frequency (single-shot) 
— Potential for time resolved velocity 
— Less accurate — estimating Doppler 
shift with few measurements 


lin 


Gon ee ea 
oo ° Jats 


Signal Level sy 


e shape function glyv,P, T) 


absorption 


ar line shape 


downshifted upshifted 


laser frequency v 
pis ab 


 ¢ laser sheet =, ae sheet Ob = 


ye ba 


un a aa 


125.66 225.662 


', “Review of 
)-Speed 
x ‘ospace 


I 
0 160 320 480 


B. Hiller and R. Hanson, “Simultaneous planar measurements of 
velocity and pressure fields in gas flows using laser-induced 
fluorescence,” Applied Optics, 27(1), January 1, 1988. 


| Doppler Shift Velocimetry Ss 


¢ Benefits 


— Capable of providing 1- or 2-component velocity maps 
— Relatively easy to implement 


¢ Drawbacks 
— Typically not time resolved 
— Typically requires reference measurement 


— Spectral shift resulting from increased pressure needs to 
be accounted for 


— Requires a gas ‘seed’ to either be naturally present or 
introduced into the flow 
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Molecular Tagging Velocimetry 


Single-Laser MTV S& 


Biacetyl (Hiller et al, 1984; Koochesfahani, 
2007) 

Acetone (Lempert, 2002) 

NO (Danehy, 2003; Hsu, 2009; Bathel, 2011) 
|, (Balla, 2013) 

N, using FLEET (Femtosecond laser electronic 


excitation tagging) (Michael, 2011) 
— Does not require seeding 
— Does require an expensive fs laser 


LIF lifetime must be long enough so molecules 
still emit light after (slightly) moved 


6) LIF Molecular Tagging Velocimetry (MTV) Ss 
Single-laser LIF MTV 


3 detect 
a 


NO naturally 
occurring in 
this flow 


— Boundary layer: slow moving gas near surface 
¢ Important: heat transfer, drag, separation, etc. 
— Instantaneous measure of 1 velocity component 


¢ Can be extended to 2 or 3 components 


P. M. Danehy, S. O’Byrne, A. F. P. Houwing, J. S. Fox, and D. R. Smith, “Flow-Tagging Velocimetry 


for Hypersonic Flows Using Fluorescence of Nitric Oxide,” AIAA Journal, 41(2), February 2003 


MTV: Delayed Images Gy 


Delay = 0 250 ns 500 ns 750 ns__In(data) 


: shown here 
¢ Measurements performed on successive tunnel runs 


P. M. Danehy, S. O’Byrne, A. F. P. Houwing, J. S. Fox, and D. R. Smith, “Flow-Tagging Velocimetry for Hypersonic Flows 
Using Fluorescence of Nitric Oxide,” AIAA Journal, 41(2), February 2003 169 


@  Single-laser MTV 


\ / 
Tunable Laser ps i 


Sheet-forming lenses Ao 


Flow Direction 


¢ - 
+ 

¢ 

s 


Molecules move, fluoresce 


Camera captures image 
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| Sin 


gle-laser MTV 


=| ZT “fF a10.0 7 
ost oat a 03 ‘ a ost 
= = | AA = 02 “A g o.2t 
5 02f & 02 Ae = i = 
a . Pe ( o.tf i o.1f 
Asi — i. ——7000° 0 
Se teil Tested Wa : Nuon) Z 
Spa 0.00cm#) [Tip yt 
| Hee - - 7 
oe ee ~ oN 
SS ino - . . 
en ee 
: 5o™ _ 1 , & 
3N) 3U1 30) 
¢ Velocity Measurement 


Instantaneous (<1 us), single-component measure of 25 profiles simultaneously 
Developed thorough uncertainty analysis methodology 


Data used in subsequent CFD study by lyer et al. of Univ. of Minnesota (AIAA 
2011-566) 


B. F. Bathel, P. M. Danehy, J. A. Inman, S. B. Jones, C. B. Ivey, C. P, and Goyne, 
“Velocity Profile Measurements in Hypersonic Flows Using Sequentially Imaged 
Fluorescence-Based Molecular Tagging,’ AlAA Journal, 49(9), September 2011 
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@ Single-laser MTV 


¢ Analysis 7 

¢ Uncertainties in mean: as good as 
3% of freestream (~1000 m/s) for 
38 images 

¢ Modeled error from intensity decay, 
exposure duration, motion blur 

¢ Error from velocity component 
parallel to MTV profiles quantified 


- 
N 


——No Decay 


eee » Decay 


Errors up to ~£5% of | easured 


0.06 


Scaled No-Decay Intensity [-] 
- 3 
rv) 
ie) 
> 
< 
©. 
ie) 
o 


Scaled Dpcay Intensity [-] 


¢ Study conclusions 

¢ Velocity gradients can bias 
measurement 

¢ Scatter off of the model surface 
limited proximity of measurement to 
the wall 

¢ Limited to low-pressure conditions 

¢ Need higher pressure capability for 
higher Reynolds number required to 
study laminar-to-turbulent transition 
at least on this short model 


d single-laser MTV: 


¢ Simple principle: 
— Works with N, seed 


— Femtosecond laser is focused down ee 
to a line or spot spas 

, -2 0 2 4 6 8 10 12 

— Gas glows for 10’s of us : = 

. 234 
° Originally developed at inp A, 
. 32 

Princeton: I 
— Air jet at room temperature, as 
atmospheric pressure oi 

— Non-linear process aaa ee 


¢ Recently applied to transonic 


f R. Burns, P. Danehy, S. Jones, B. Halls, and N. Jiang, 
cryogenic flow at NASA Lan gley “Application of KUBET Velocimetty in,thle NASA ney Oe ce 
Meter Transonic Cryogenic ICO 
. ‘a 31SLAIAA. Aerodynamic I 
— Tracking of spot provided 2-D Ground Taste 


velocimetry 


— ~1 m/s accuracy and precision 
J. Michael, M. Edwards, A. Dogariu, and R. Miles, “Femtosecond 
laser electronic excitation tagging (FLEET) for quantitative velocity 


imaging in air,” Applied Optics, 50(26), September 10, 2011. Me 


| Single-laser MTV S& 


¢ Benefits 
— Instantaneous measurement 
— Capable of providing velocity information along a line (1- 
component) or at grid points (2-component) 
— ‘Seedless’ methods available 
— Relatively simple to setup 


¢ Drawbacks 


— Reduced spatial resolution compared to Doppler shift 
approach 

— Operational limits often dictated by fluorescence lifetime 

— Potential for model damage 

— Potential for altering thermodynamic state of gas 
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6 Multilaser MTV S& 


¢ Typically, one laser generates a line, series of lines or 
a grid pattern by, ionization (LEI*2°72°), vibrational 
excitation (RELIEF*¢? +3"), photolysis, bleaching or 
other process. Additional laser(s) read out pattern. 


— Example: Excimer laser at 193 nm+N,+O,9NO 


— Air Photolysis and Recombination Tracking (APART) 
¢ Dam et al, Optics Letters, 26 (1), pp. 36-38 2001 


¢ A list of partner species used in the photolysis writing 
and reading process includes: 
H20-OH, 33-139 NxO-NO, 4° O2-O3, 41-142 
N2/O2-NO, !43-146 and NO2-NO. 116-147-154 
¢ Method can use air constituents (N,/O,, O,, N,) or 
can be seeded (H,O, N,O, NO,, Na, etc.). 
— We have used NO,NO to study transition 
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@ Multiiaser MTV 
> NQd oa khetabrris 


Probe Beam 


a 


— == Pump Beam 


(for details of velocimetry method, see Bathel et al, AIAA-2011-3246 and others incl. TAMU'*2, OSU"*°, etc.) We 


0 250 500 750 1000 1250 1500 
Streamwise Velocity (m/s) 


0 250 500 750 1000 1250 1500 
Streamwise Velocity (m/s) 


¢ Transition-to-turbulence study (Bathel et al. AIAA-2011-3246) 
— Reynolds number requirement necessitated higher operating pressures 


— Signal lifetime too short for single-laser NO LIF MTV 
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| 3 laser MTV: measure u, v, T S 


¢ Texas A&M has used 3 laser technique with crossed beam pattern 
for 2-component velocity and temperature in an underexpanded 


N, jet flow with NO,: 
— VENOM: vibrationally excited nitric oxide monitoring 
— Excite to v =1 state 


400 ns delay 800 ns delay ‘ measured 


Raw data images U-Velocity (m/s) Temperature (m/s) 
(~5% uncertainty) (~9-35% uncertainty) 


R. Sanchez-Gonzalez, R. Srinivasan, R. D. W. Bowersox, and S. W. North, “Simultaneous velocity and 
temperature measurements in gaseous flowfields using the vibrationally excited nitric oxide monitoring 
technique: a comprehensive study,” Applied Optics, 51(9), pp. 1216-1228, March 20, 2012 
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d Multilaser MTV: Unseeded S& 


Air Photolysis And Recombination nae 
Tracking (APART) 
— Nitric oxide generated in unseeded 
air via dissociation by UV excimer 
laser beam Tag 
— NO then probed with LIF excitation 532 and 580 nm 
wavelengths 


¢ Raman Excitation and Laser-Induced 
Electronic Fluorescence (RELIEF) 
— O, vibrational excitation via Raman 
pumping 
— UV Excimer laser used to probe 1** 
vibration state of oxygen 
— Works for flows where T < ~750 K 


coy i a ee 


mixing aed Ata ek 


Dye Cuvette 
to Seed 


Raman 


tat Epa ea 


ArF Excimer Cylindrical Lens 


Laser 0 
193 nm 


Cylindrical Lens 


Test Cell 


CccD 
ND:YAG LASER 
(SECOND CAMERA 


HARMONIC) 
32nm 


2 > Focusing Lens 


1a 
Oxygen Raman Cell 
1st Vibrational Q Branch 
532 --> 580nm 


Cell 


é Multi-laser MTV S& 


¢ Benefits 


— Capable of providing instantaneous velocity along a 
line (1-component) or at grid points (2-component) 

— Several ‘seedless’ methods available 

— Not limited by fluorescence lifetime (works in high P) 


¢ Drawbacks 


— Reduced spatial resolution compared to Doppler shift 
approach 


— More complicated setup — requires multiple lasers 
— Potential for model damage 


— Altering thermodynamic state of gas 
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Summary 


| Velocimetry Techniques Ss 


¢ Introduced three fluorescence-based techniques 
— Doppler shift velocimetry 
— Single-laser molecular tagging velocimetry 
— Multi-laser molecular tagging velocimetry 


¢ Provided representative examples for each of the 
techniques 


¢ Discussed benefits and drawbacks for each approach 


¢ Following presentation will provide a more in-depth 
discussion of laser-induced fluorescence used for many of 
these MTV techniques. 
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Backup Charts 
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Experimental Setup: Top View 


T2 Free Piston Shock Tunnel, ANU 


Thermometry laser sheet 
( 5 mm wide, 0.7 mm thick) 


Velocimetry laser sheet 
(3 mm wide, 0.2 mm thick) 


Velocity ~ 3000 m/s; Mach Number = 7, laminar boundary layer 
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o Average of Single-Shot Measurements S 


Height Above Flat Plate (mm) 


16 


—— Average Velocity 
12 Standard Deviation 


10 


—_ 
a 


6 Uno = ~3000 + 100 7 


0 ————<——__ 
0 500 1000 1500 2000 2500 3000 3500 
Velocity (m/s) 


Freestream velocity is uniform at 3035+60 m/s. 
95% boundary layer thickness is 1.4+0.1 mm. 
Single shot measurement precision: ~60 m/s out of ~3000 m/s (2%) 
- spatially averaged uncertainty (freestream region) ~100 m/s (8%) 
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| Three-Laser MTV (NO,~>NO velocimetry) 


(for details of velocimetry method, see Bathel et al, AIAA-2011-3246 and others incl. TAMU'49, OSU'*®, etc.) 


—_— | PLIF 


me Cart B 
PLIF Cart A 


NO,+355 nm > NO 


Write lines of NO into 
NO, seeded in the flow 


6 Flow Visualization: No trip 


¢ Flow appears laminar when no trip is present 


| Velocimetry: 1-mm Tall Trip (centerline) BY 


Waill-Normal Pesition (rar) 


5 


K 


& 


Streamwise Posi ion (mm) 
Streamwise Position (mm 


Velocity profile is full near flat plate with a deficit near edge of BL 


Note: 2x lower flowrate in this run compared to no-trip case previously shown. 
U’ data has increased to peak ~250 m/s (3x larger than laminar case). 


Shape of U’ profile changes slightly with distance downstream. 
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Velocimetry: No Trip 


flat plate 
24 
2B 125 m/s : 3 
2.9 <10005/s 3 ! 
18 3 3 
16 | | ! 
Bis ; S ? 
£12 i : i 
1.0 | : 
0.8 oN 4 1 ny i 2 1 
y ip i 
£06 ww gh st j wv i 4 
04 i =A a } 4 : y r* Pa 
: I 4 X ry ' i : 
02 yee rit re Fon ar Zi 
0.0 LA i Me < 1 1 a, v 
eS ee secsesssssssseePReBSeseese ess SERRE 
ad we Gy os Be hao we Sb ois Wins ae ae os WwW © WH &w B OW 


¢ Data formatting: SERINE REISE RY 
— Center of the data point is mean velocity; Width is uncertainty in mean 
— Standard deviations shown in red (instrument error ~2% of max U was subtracted) 
¢ Blasius-like boundary layer profiles observed (to within 100 um of wall) 
— Deficit near edge is likely caused by seeding in zero velocity gas into BL (see Johansen et al) 


— U’~80m/s at maximum. Peak U’ is ~ 0.7 mm above surface. 
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1-mm Tall Tr 


| Velocimetry 


"113.8 


Streamwise Position (mm) 


4.5 mm profiles show profiles most different from laminar 


Flow becomes more laminar far from centerline (similar to flow vis) 


U’ drops back to ~80 m/s 
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measurements 
on centerline 
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int. 
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Streamwise Position (mm) 


0.53 mm trip profile is between no and 1-mm but no inflect 


measurements 
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| Comparison of U’ for different trip heights 
No Tr 
0.53-mm Trip 


Streamwise Position (mm) 


1.0-mm Trip 
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+ 114.3 
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108.6 
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¢ 0.53 mm trip U’ peak much closer to plate. 
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é Three-laser MTV: uncertainty & design S 


¢ Analysis a 
¢ Uncertainty for high s/n data 
- Mean: 1% of edge velocity 21s 
Single sho 1.2% of edge velocty KA 
- Single-shot: 1-2% of edge velocity 


¢ Measurement within 100 um from wall 


” @ precision 
- Tradeoffs for experiment design = oe? fbesttypical 
- Accuracy and precision improve with = § owe 
probe delay 3 1m u 
- Spatial resolution worse with longer = et? 
probe delay > rf) @ 
- Can use longer delay for slow flow  & is a 
* Current configuration requires < | oO » |° 
relatively large blowing rates for : f / == 7 
adequate signal-to-noise 2 adeutaby Lis = 
- May be significant wall-normal S$ (best/typical) 
velocity component 
1 10 100 


Probe Laser Delay, Atproge (US) 
Data obtained at static, low pressure conditions !% 


| Intro to Rayleigh and Raman Spectroscopy for in seal 
Hypersonic Nonequilibrium Flows 


Paul Danehy, NASA Langley Research Center, Virginia, USA 
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é Rayleigh and Raman Scattering S& 


a 


a 


¢ When light passes through a gas: 
— LIF is a resonant absorption/emission process 
— Mie, Rayleigh and Raman are scattering processes 
¢ Mie Scattering: d=A (dis the particle size) 
¢ Spontaneous Rayleigh and Raman scattering: d << A 
¢ Will explain Rayleigh and Raman two ways 
— Physical interpretation 
— Classical description 
— (Not going to discuss quantum treatment) 196 


é Rayleigh Scattering Ss 


¢ Rayleigh scattering: 
— Electric field causes electron cloud to oscillate 
with respect to the nucleus: induced dipole 


Vv.—~>e@e2@e@e8- Vv, 


http://www.timkelf.com/Research/ResearchSERS.html 


— An oscillating dipole acts like an antenna: radiates! 
¢ Radiation is perpendicular to oscillation 
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Many Particle Scattering Ss 


¢ A plane-wave perturbs many particles: they radiate 
— Interfere constructively in direction of propagation 


— Interfere destructively in other directions 


_ NO SCATTERING 


E ~ 
INCIDENT OF on ay ; 
Local clear zatlon =e 


WAVE 
Sources \ 


eet f wae a we COHERENTLY 
ail iw i | “ rs eel SCATTERED 
i | (i as | 
/ / x { FORWARD 
| > 
/ \ } a RADIATION 
/ \ 
/ \V 
/ ea 
H j = ie 
| 
\ ; 
\ |! G.C. Baldwin "An 
\ iN. WAVELETS AT 90° introduction to nonlinear 
optics." New York: 


INTERFERE DESTRUCTIVELY 


Plenum Press, (1969). 


— |f exact same number of particles everywhere in gas, 


destructive interference would give no side scatter 
¢ But random statistical variation in particle concentration: 


deconstructive interference is incomplete: Rayleigh Scattering! 198 


é Rayleigh Scattering S& 


¢ Very weak process 
— 1 photon in 0.1 million scatters / meter at STP 


¢ Wavelength dependence: 1/A7 


Rayleigh scattering in the atmosphere after sunset, 
picture taken over the ocean, at 500m altitude. 
http://ms.wikipedia.org/wiki/Selerakan_ Rayleigh 


¢ Causes blue sky and red sunsets 
— Mie scattering from aerosols also contributes 
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| Raman Scattering Ss 


¢ Similar to Rayleigh scattering but molecules gain or 
lose energy during scattering process 
¢ Raman x1000 lower probability than Rayleigh 


virtual 
energy Energy Level Diagrams 


level 


energy molecule 


stays same gains energy molecule 
loses energy 


tA bb) 8 


N, molecule 


o Raman Scattering S& 


¢ Raman scattering results in a frequency shift 
— Molecule vibrating at frequency v 


mol 


Rayleigh Scattering 
V. > @e@e@e8-— \, 
Raman Scattering 


ad 1 hee 
ed 


V Modulates the frequency 
mol 


http://www.timkelf.com/Research/ResearchSERS.html 


201 


| Rayleigh and Raman Scattering S& 


Oy a ee es ee ee es Se ee ee 
Po. lad ee hy RE a tay 
I, , Stokes 
Rayleigh aaa, anti-Stokes 


Raman 


¢ Energy level diagram: 
— Stokes Raman: energy absorbed by molecule 
— anti-Stokes Raman: energy lost by molecule 
¢ Initial state must not have been a ground state 
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| Rayleigh and Raman Spectra Ss 


Rayleigh 


\ Stokes 
anti-Stokes esapaal vibrational 
vibrational por eneene Raman 


Raman 


Maser Wavelength, A 


¢ Rayleigh: near laser’s wavelength 
¢ Raman: frequency shift depends on molecule, 
State (sensitive to concentration and temperature) 


— Rotational Raman 
— Vibrational Raman 
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é Rayleigh and Raman Theory 1 Ss 


¢ Scattered intensity, /, is proportional to the 
square of the induced dipole moment, P : 


B= ce ,ak 
where &j is the permittivity of free space, o is 


the molecular polarizability and £ is the 
incident electric field given by: 


E = E, cos(wy t). 
where E, is the amplitude of the electric field, 


@, is the frequency of the laser light and f is 
time. 204 


| Rayleigh and Raman Theory 2 Ss 


¢ Polarizability, a, depends on internal structure 
— varies with time during vibrational oscillations at 
the natural frequency of the molecule, @, 


O—O —— 0 


— a varies with spatial coordinate, O 
¢ More loosely held electrons > more polarizable 
¢ Taylor Series Expansion: 


0 
Oo = Ag + (53) dQ 
0 


where dQ is small displacement about equilibrium 
positions which occurs during normal vibration 


dQ = Q,)cos(w,t) 
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6 Rayleigh and Raman Theory 3 S& 


¢ Substituting previous equations into p = SOE: 
a= |. + (<=) Qo cos(wyt)} SoEo COS( Wot) 
AQ/ 


Expanding and using trig identity: 


- —> O — 
DP = ApEgE_ COS(Wot) + (53) EQ Ep|COS(W— — Wy)t + COS(W— + W,)t] 
0 


Rayleigh Scattering F— ——~=~~~~Stokes Raman [— anti-Stokes Raman 
Resulting induced dipole moment (and 
scattering) is at three frequencies 


¢ Amplitude is usually determined empirically: 


(2) _e 
dQ/,, NI 


Differential Cross Section scattered power irradiance 


(varies by molecule and with | but not with P, T) per solid angle diff. cross number 
(zz indicates polarization of laser and detection) section density 206 


| Rayleigh and Raman Theory 4 @y 


¢ Selection Rules: Raman Activity 
— Not every molecule or vibrational mode produces 
Raman 


= — 0 — 
D = ApEgEy COS( Wot) EG) or [cos(Wo — W,)t + COS(Wy + W,)t| 
4 


— Amplitude of Raman is proportional to the 
magnitude of: ~< 

— Polarizability, «, must change during vibrations to 
be Raman active 


— Raman active: —_ 


N 


diatomic \ linear triatomic 
symmetric 
stretch mode 


— Raman inactive: 


linear triatomic 
Ems - 
asymmetric 
stretch mode 207 


| Rayleigh and Raman Theory 5 S& 


¢ Rayleigh Scattering is at or near the laser’s 


wavelength. 
— Doppler broadened (~T’”) 


— Can be Doppler shifted (~V) a 
fe = ran NI 


— Signal proportional to N (~o) aO 
Shift, V 

| 
=) | 
= 
F Amplitude, p 
£| | 

Maser Wavelength, Xd 


— Can potentially measure these simultaneously 
¢ Requires high spectral resolution 208 


| Rayleigh and Raman Examples Ss 


¢ Rayleigh 
— Velocimetry, multi point 
—V,T, p, single point 
—V, T, p, imaging 
— 9 measurement (and infer T) 


e Raman 
— Concentration 
— Temperature 
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| Rayleigh Velocimetry Ss 


¢ Rayleigh scattering velocimetry can measure 


selected velocity components 
— Measured component is bisector between 
incident laser and collection optics 


k Wave vector: 
seme k=27/ 


D. Bivolaru, P. M. Danehy, R. L. Gaffney, Jr. and 
A. D. Cutler, “Direct-View Multi-Point Two- 
Component Interferometric Rayleigh Scattering 
Velocimeter,” 46th AIAA Aerospace Sciences 
Meeting and Exhibit, Reno, Nevada, AIAA Paper 
2008-236, January 7-10, 2008 


— Can measure multiple velocity components by 
adding additional detection angles or laser beams 
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| Rayleigh Velocimetry Ss 


¢ Rayleigh scattering velocimetry has been 
performed different ways, always requiring a 
high resolution device to resolve the Doppler 


shift 
— I, filter in front of CCD | 
¢ Imaging capability | 
¢ limited dynamic range . 
— Fabry Perot etalon = aemoane 
¢ PMT based (faster) 
¢ CCD based (multi-point) 
— Combined J, filter and Fabry Perot etalon 


Shift, V 


Amplitude, p 


Intensity 


— See manuscript for references 211 


Rayleigh Velocimetry using Fabry-Perot Etalon S& 


¢ Velocimetry using a Fabry Perot Interferometer 
(FPI) with CCD camera 


collimated 
input beams 
L, etalon Mae 
~ PBS Lg M 
a Ls 


em CCD 


Ly 
OF 


a @ 532 nm 


measurement 
volume 


reference Rayleigh signal 


° D. Bivolaru, P. M. Danehy, R. L. Gaffney, Jr. and A. D. Cutler, “Direct-View Multi-Point Two-Component 
Interferometric Rayleigh Scattering Velocimeter,” 46th AIAA Aerospace Sciences Meeting and Exhibit, Reno, 
Nevada, AIAA Paper 2008-236, January 7-10, 2008 212 


Rayleigh Velocimetry using Fabry-Perot Etalon 


one-component system 


measurement 
volume __ 


wa 
- —= 
a 
- 


= reference = 
a 


collection = 
F-P Etalon is two = 
planar, separated 


optics ww 
x 
. ~~ mirrored surfaces 
etalon 
> 


Spatial and Spectral Info 
~ > 
<i 


bd 


Tk en an 


D. Bivolaru, P. M. Danehy, R. L. Gaffney, Jr. and A. D. Cutler, 


“Reference Bivolai Dar ( 
Wg ana “Direct-View Multi-Point Two-Component Interferometric 
Frequency Rayleigh Scattering Velocimeter,” 46th AIAA Aerospace Sciences 


Meeting and Exhibit, Reno, Nevada, AIAA Paper 2008-236, 
January 7-10, 2008 213 


Rayleigh Scattering Image Processing Ss 


Raw Image data D. Bivolaru, P. M. Danehy, R. L. Gaffney, Jr. and A. D. Cutler, “Direct- 
View Multi-Point Two-Component Interferometric Rayleigh Scattering 
Velocimeter,” 46th AIAA Aerospace Sciences Meeting and Exhibit, 
Reno, Nevada, AIAA Paper 2008-236, January 7-10, 2008 


Import to Matlab® 


x-y to r-y transform 


Circles provide reference 
(7.5 GHz etalon) 


Combine to spectrum and fit 


¢ Significant image 
processing is required 
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@) Rayleigh Scattering Velocity Results €3Y 


Normalized Intensity 


0-44 0.5 4 45 1.9 
Radius, a.u. 


¢ Rayleigh scattering in a large scale (6 cm diameter) Mach 1.6 H,-air 
combustion heated jet operating at Mach 5.5 enthalpy at NASA Langley 


Velocity precision of ~40 m/s out of 1200 m/s, (3%) 
The dynamic range of the instrument was ~3000 m/s 
- precision 1% of dynamic range 


D. Bivolaru, P. M. Danehy, R. L. Gaffney, Jr. and A. D. Cutler, “Direct-View Multi-Point Two- 
215 Component Interferometric Rayleigh Scattering Velocimeter,” 46th AIAA Aerospace Sciences, , 7 
Meeting and Exhibit, Reno, Nevada, AIAA Paper 2008-236, January 7-10, 2008 


Rayleigh Scattering: T, V, 0 Ss 


Image 
Dissector 


Fiber 
Image 
Edge 


PMT 3 PMT 4 


¢ Continuous laser excitation 
— 10Wat532 nm 


e Fiber optic collection 
¢ Data rate: 10 kHz 
e Single point measurement 


A. F. Mielke and K. A. Elam. "Dynamic measurement of temperature, velocity, and density in hot 
jets using Rayleigh scattering." Experiments in fluids 47, no. 4-5 (2009): 673-688. 2G 


Rayleigh Scattering: T, V, 0 S& 


a 1.50 
0.55mm dia, 75-m | 
length optical fiber 1.254 
to spectral analysis | 
and detection system 1.004 
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~ 
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near-sonic jet at NASA —* 


0.254 


Glenn Research Center A A__| SA eer —__| 


x/D x/D 
A. F. Mielke and K. A. Elam. "Dynamic measurement of temperature, velocity, and density in hot 
jets using Rayleigh scattering." Experiments in fluids 47, no. 4-5 (2009): 673-688. 
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| Rayleigh Velocimetry (J, filter) Ss 


¢ Instead of an etalon, use I, filter > 2D imaging 


1.00 pA 


|W ilo 


|, absorption 
near 532 nm 


Pressure 
broadened 
absorption 


Background 


normalized transmission 


: filter 
scattering * ‘A 

0.40 | a 

f 

} t 

/\ Doppler 

0.20 = shifted 

[ scattering 

1 
0 00 i... it W | 


Vv 
G. S. Elliott, M. Samimy, S. A. Arnette “A molecular filter 
18787.5 18788.0 18788.5 based velocimetry technique for high speed flows,” 
R. B. Miles, W. R. Lempert and J. N. Forkey, ay Experiments in Fluids, Volume 18, Issue 1-2, pp 107-118, 
“Laser Rayleigh Scattering,” Meas. Sci. Technol. frequency (crm) December 1994 
p. R33-R51, v.12, 2001 


° Gas cell, filled with low-pressure N, plus some I, 
crystals; control temperature to determine absorption 
¢ Scan Nd:YAG laser frequency (or fixed) 218 


@ Rayleigh Imaging of 7 V; o(or P) S& 


Streamwise distance 


from nozzle exit (mm) 


¢ Direct YAG laser sheet into a Mach 2 jet 


j 170K >, 1050 torr ’ 235 m/s 
23 S = a s e are! 
gu “ gwar 22 
“ux az 
“= — >» ie if ae 
“UP oe : Bi a . 
2 160K ao 2 850 torr a2 3 220 m/s 
Ea = 8 
ag aS 
g = a — 
e¢« g« 
na 2 n £ 
150K 650 torr 205 m/s 
32h? 2 3 i; oT 2 
Spanwise distance Spanwise distance Spanwise distance 
from jet axis (mm) from jet axis (mm) from jet axis (mm) 


R. B. Miles, W. R. Lempert and J. N. Forkey, “Laser Rayleigh Scattering,” Meas. Sci. Technol. p. R33-R51, v.12, 2001 


¢ Scanned laser frequency, determine 7 VY. o (report P) 
based on measured spectrum at each ah 
¢ Same method can be used for velocimetry from 


particles (usually fixed frequency, multiple detection angles) 
— Planar Doppler Velocimetry (PDV) or Doppler Global Velocimetry (DGV) 
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Rayleigh Scattering: 0 Ss 


Shift, V 
Rayleigh scattering has a 7 | 
spectrum: cl | Amplitude, p 
— Sensitive to 7, V; p =| | 
— Prior examples used Fabry Perot Mieser Wavelength, A 


Etalon or I, cell to observe spectrum : 
— These devices reject signal, making 
them less sensitive detectors of ¢. 


Collect al// of the scattering 
— Less equipment/less setup/simpler ° iobeaipeiceees 


rror bars 


a 


e 
oO 


_ 
a 


Rayleigh Density (x1017/cm3) 


= 
o 


— Bigger signal ” 0.5 mantels Daivsiiy 0am 4.0 4.5 

— Cando single point or imaging Using a UV laser (e.g. ArF at 193 
ae nm) allows measurements to be 

= Only sensitive to Pp: made at low pressures (e.g. 1/300 


wy atmospheric density) 
O (If you know P and composition you Balla, R. Jeffrey, and Joel L. Everhart. 


can convert p to 7) "Rayleigh scattering density 


measurements, cluster theory, and 
nucleation calculations at Mach 10." 
AIAA Journal 50, no. 3 (2012): 698-707-~" 


é Raman Spectroscopy asa 


¢ Raman spectrum obtained in a CH,-air flame: 


— Average of 500 pulses using stretched 0.5 J/pulse, 532 nm 
¢ Stretch 10 ns pulse to 500 ns to avoid laser-induced breakdown 


2 N, 
3 6000 
~. 
> Pure-rotational 
© 4000 . — 
e anti-Stokes ands 
re N 
£ 2000 2 
0 
450 500 550 600 650 700 


Wavelength (nm) 


Figure courtesy of Jun Kojima, NASA Glenn Research Center (Jun.J.Kojima@nasa.gov) 


Can measure species (N,, CO, CO,, H5, HO) from 
such spectra &temperature from Stokes/anti-Stokes 
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Raman Concentration Measurement 


ww _(and Rayleigh Temperature and LIF concentration) 


¢ Example of Raman concentration measurements in a subsonic 


turbulent CH,/air flame at Sandia National Laboratories, California 
— 4Nd:YAG lasers at 532 nm to produce 1.8 Joules/pulse, focused to a point 


0.08 


Rayleigh . 
camera =< — 


300 600 900 1200 1500 1800 2100 
T (K) 


R. S. Barlow, G.-H. Wang, P. 
Anselmo-Filho, M. S. Sweeney, 
and S. Hochgreb, “Application 
of Raman/Rayleigh/LIF 
diagnostics in turbulent stratified 
flames,” Proceedings of the 


0.60 0.60 Combustion Institute, 32, pp. 
300 600 900 oe 1500 1800 2100 300 600 900 > 1500 1800 2100 945-953, 2009 


x; precision ~0.7% to 7.5% 
T precision ~0.75% 


7 Raman/Rayleigh/LIF line imaging (~10 pixels / mm along the 6-mm-long probe volume) 
N,, O05, CH,, CO, H,0 and H, from Raman, CO from 2 photon LIF 


— Temperature from Rayleigh density measurement (idea gas law) eee 


| Raman Temperature Measurement Ss 


¢ Several ways to measure temperature with Raman 
— none as precise as CARS (3%) 


Rayleigh 


Stokes 
anti-Stokes aieaal vibrational 
vibrational nhoLion Raman 


Raman 


Raman 


Maser Wavelength, A 


¢ Rotational Temperature from pure-rotational Raman 


¢ Vibrational temperature measurements eee nnavenwe 
— Only good above 700 K 
— Stokes / anti-Stokes vibrational Raman 
¢ Low resolution, so can measure concentrations at same time 
¢ ~20% precision in measurements along a line (Wehrmeyer, 1996) 
— Spectrally resolved vibrational band (Sharma et al, 1993) 
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Raman Temperature Measurement Ss 


¢ Raman temperature measurement from spectrally 


resolved vibrational N, band. 
— Excited with KrF Excimer laser 0.25 J at 248 nm (25 ns pulse) 


T,=3500 K 


2 
e 


= 
N 


. on 
@ = T, Experiment “Sine: 


Photon Counts 
ad 
oa 


Normalized Temperature T/T, 


ad 
a 


262.2 262.4 262.6 262.8 263.0 263.2 263.4 263.6 2 0 2 
Wavelength, nm : Nozzie Position, cm 


_ 4 6 


S. P. Sharma, S. M. Ruffin, W. D. Gillespie and S. A. Meyer, “Vibrational Relaxation Measurements in an Expanding 
Flow Using Spontaneous Raman Scattering,” J. of Thermophys. and Heat Trans. Vol. 7, No. 4, Oct.-Dec. 1993 


— Single shot precision ~5% (not explicitly stated) 


¢ More precise than Stokes/anti-Stokes but cannot measure species 


concentrations with same instrument. ae 


C) 


Rayleigh Scattering: 0 S 
¢ When precise measurements are required, 
high S/N is needed 


— Use high powered lasers, focused to a point 
° If too high: laser-induced breakdown 


¢ Sandia, Livermore, Calif.: 


Raman/Rayleigh/LIF line imaging at P = 1 


4 double-pulse Nd:YAG lasers: 1.8 J/pulse * 


Camera images Rayleigh along line 


atm 


— ~10 pix/mm along 6-mm-long probe volume 
— N,, O,, CH,, CO,, H,O0 and H, from Raman 


CO from 2 photon LIF 

Temperature from Rayleigh density, p, 
— Need composition to determine ( 2% 
— If known P, idea gas law give T \@Q 
— T precision reported: ~0.75% 
— (pprecision is probably same) 


~~) 2 
0 


Turbulent CH,/air flame 


R. S. Barlow, G.-H. Wang, P. Anselmo- 
Filho, M. S. Sweeney, and S. Hochgreb, 
“Application of Raman/Rayleigh/LIF 
0.65 diagnostics in turbulent stratified flames,” 
Proceedings of the Combustion Institute, 
32, pp. 945-953, 2009 


0.60 
300 600 900 1200 1500 1800 2100 
T (K) 


Cutting edge MHz laser: 40x lower energy (Jiang, 2011) 


é Advantages and Disadvantages: Rayleigh & Raman S& 


¢ Advantages: 
— No seeding required 
— Only one color of laser used (simplicity) 
— Multiple properties measured at once 
° TV, efor Rayleigh 
e T, many species for Raman 
¢ Can do Raman and Rayleigh together and with LIF 
¢ Disadvantages: 
— Weak signal: 


¢ Requires very high power laser 

— stretched to avoid laser induced breakdown 
¢ Requires large windows, big, close collection optics 
¢ Usually singe point or line measurements 


— Temperature less accurate & precise than CARS 
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Intro to Coherent anti-Stokes Raman Spectroscopy 
(CARS) for Hypersonic Nonequilibrium Flows 


Andrew Cutler, The George Washington University, USA 
Paul Danehy, NASA Langley Research Center, Virginia, USA 
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Outline 


Brief introduction to nonlinear optics 

Introduction to CARS and basic theory 

NASA / GWU CARS System 

— Temperature and composition measured in UVa 
dual mode scramjet 

— Thermometry shows vibrational nonequilibrium 

CARS in a non-thermal plasma 

Concluding remarks 
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| Nonlinear Optics 


First, a few words about nonlinear optics 
— Conservation of energy 


— Conservation of momentum 
— Wave equation from Maxwell’s equations 


¢ Second Harmonic Generation as an Example 
— 2° order nonlinearity (solid) 


— CARS is a 3 order nonlinearity (gas) 
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| Intro: Nonlinear Optics Bes 


— Second Harmonic Generation (SHG) 


1064 iN 
nm 1064nm 532 nm 


— Solved using the Wave Eqn: 


— Assume SVEA and 1° Born Approx: 
D0 4 
a PAL 
> A? 
—disacoefficient, P is power, L is length, A is area 
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G.C. Baldwin "An 
introduction to nonlinear 
optics." New York: 
Plenum Press, (1969). 


Input E at w 


t 
¢ Non-linear effects: Intense electric fields of one 
frequency can generate a P at another 
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| Phase Matching 


¢ What is Phase Matching? 
— Ensuring that the light wavefronts ‘match up’. 


__———_ Fundamental 
Phase nj [\— Py vir harmonic 
matched 
Not phase DASA LAY |k| = OTA 


matched =on/Cy 


° ven, Causes Sen interference 
¢ New photons generated are out of phase with previous 
ones: resulting in loss 
¢ For SHG there are tricks to make it phase matched 
¢ Conservation of momentum: phase matching 


¢ Conservation of energy: ®@,=20, —=——=——Ci‘é: 
2 1 


o What is CARS? 


Broadband N, CARS 


- 

— 
oe 

=—— 


Interaction 
Region 


Spectrally-broad red beam + 
Spectrally-narrow green beams = 
single-shot blue CARS spectrum 


@ CARS Energy Level Diagram Gw 
w Ss (conservation of energy) 


e Lasers: 2 Green (532 nm) and 1 Red (607 nm) 
— Frequency difference between Green and Red 
Probes N, vibrational Raman Shift 


iS 


¢ Result: Blue signal beam at 473 nm. 


— Spatially, spectrally, temporally filter this signal 
¢ Good technique for luminous flows 


— Send to spectrometer equipped with CCD camera 
— Measured CARS spectra is best fit to a library of 
theoretical soectra on a single-shot basis (10 ns) 
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CARS Phase Matching Cw 
ww __(conservation of momentum) SHES 


ee re ee ae Planar Boxcars Geometry 


; em 
anti-Stokes 
<a _a™ 


(blue) 


N, first vibrational state 607 n 473 nm 
Ground state _ 


It is straight forward to compute the CARS phase matching angles: 


N2 cars adjust inpt 

Different molecules 
require different angles 
(can get multiple signals 
at different spatial 
0023576039 —~«|OCations: Dual-Stokes) 


N2Raman= 2322.528 cm-1 0.023576939 
(nominal N2 = 2330 cm-1) 3.00927E-20 


lambda wavenum angle (degrees 
stokes 607 16474.46 
pump 2 532 18796.99 
pump 1 532 18796.99 
anti stokes 473.4956 21119.52 
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é CARS Phase Matching Geometries 
k +k, =k, +k,., 


532 nm 932nm 


Collinear CARS: 


— Poor spatial resolution 
607 nm 473nm 


ae “<r 
— Oh _Q ™ 
607 n 473 nm 


(Planar) BoxCARS: 


¢ Best spatial esention (1.5 mm x 0.1 mm x 0.1 mm) but lower signals and 
affected by beam steering from turbulence 
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| CARS Theory 1 


¢ Non-linear effects: : 
— Intense electric fields of one frequency 
can generate a P at another 


= P=e (verge ty E+...) 


— 7%) =2nd order susceptibility 


(ie, no gases, liquids) 


— 7°) = 3rd order susceptibility 


¢ All dielectric media 


— Generally (yE >> y¥) E? >> 7) EB ) 


G.C. Baldwin "An introduction to 
nonlinear optics." New York: 
Plenum Press, (1969). 
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6) CARS Theory 2 


¢ Polarization generated at CARS signal frequency: 
P 8) (Wsignat) = €0X carsE (Wyump)E (Wstoxes)E (W probe) 
¢ Substitute into the wave equation. Result: 
| signal x | pump Istoxes! probe Ix CARS| [ir 
¢ To obtain large signals, want long L, high laser 
power (tight i and strong resonance 


XCARS = yr iT; + Xnr where Kj « NA; (= a), 


* [; is damping coefficient, NV is number density, 


A, is the population difference 
— N> concentration; A;> Temperature sensitivity 
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| CARS Theory 3 


¢ Instructive to graph example N, spectra: 


100 


— susceptibility 


80 — signal 


Different 
rotational 
lines 


al) 


square root(CARS sign 


qua 
Ww 
oO 


anna 
ea 


¢ N, is commonly probed with CARS because it is often 
present and spectroscopy is well known 


¢ CARS spectrum at 1 atm., 1400 K, v” = 0 band. 
— N, Q-branch (Av=+1, AJ=0) 
— CARS signal spectrum is convolved with lasers used a 


é CARS Theory 4 Ad 


¢ Effects of temperature (left) and composition (right): 


1000 


= a 
o c 
a bo 
@ ~=©10 a 
iH < 
< YY 
= 3 
8 8 
S 2 2250 2275 2300 2325 2375 3 
toy a 
a a to 79 
v’=2 v’=1 v’=0 
0.1 
Gai 2225 2250 2275 2300 2335 2350 9375 


Frequency wavenumber, cm+ 


Frequency wavenumber, cm 


Vibrational “hot band” becomes CARS ena aves 
populated above 700 K See ge eee 

the gas composition through 
resonant and nonresonant parts 
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G@ wuttispecies Measurement: Dual-Pump CARS GW S& 


Interaction 
Region 


Interaction 
Region 
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@) Dual-Pump CARS: N,, O,, H, 


* Green and Red Probe N, Raman Shift 
¢ Red and Yellow Probe O, Raman Shift 


2 
¢ Result: Two different spectral regions of 


Raman Shift appear coincident in the 490 nm 
CARS beam. 


— ‘self calibrating’ compared to other methods 


— H, pure rotational Raman lines present 
O’Byrne et al AIAA Journal v. 45, 2007 
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é Typical Dual-Pump CARS Spectra and Fits aw 


is 4 Freestream 


¢ Single-shot spectra 
obtained in combustor 


shown with best fit 

¢« Absolute measurement 
of rotational, vibrational 
temperature, N., O,, 
and H, mole fractions 


— 
[o>] 


— Experiment 
— Theory 


_ — 
yo > 


sqrt(normalized intensity) 
3 


2225 2275 2325 2375 
Wavenumber, cm* 


In fuel reacting plume 


— Experiment 
— Theory 


N 
on 


¢ Measurement obtained: 
— In UVa Scramjet (Conf C) 
— p=0.18, Plane 2, 2/H=0 
— in10 ns (instantaneous) 
— at 20 Hz repetition rate 


; — Measurement volume = 
2075 2125 2175 2225 2275 2325 2375 1.7 x 0.2 x 0.2 mm? 


Wavenumber, cm 


= 
on 


= 
oO 


sqrt(normalized intensity) 
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Set up in basement below DCSCTF 


6 Features of NASA/GWU Mobile CARS System pease 


¢ YAG laser plus 2 dye lasers on single cart 
¢ Transmission and receiving optical stages 
¢ Can put system ona truck and moved to facilities 


¢ Measures T,.., T,;,, Nz, O., H, mole fractions at 20 Hz. 
— Quantified accuracy, precision 


¢ Fully automated data acquisition 

— Remote operation of system 

— X-Y (-Z) probe volume traverse capability 

— Input test matrix; then system acquires and automatically saves data for hours 
until complete 

— Laser power automatically adjusted to prevent saturation effects (detector 
saturation, Stark broadening, Raman pumping) 

— Beam alignment monitors; remote control alignment 


¢ Beam shaping capability to maintain alignment 
measurement in turbulent/vibration flows 
— 10x improvement compared to USED CARS 
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Setup at UVa: Optical System 1 (Say 


Laserlab <—————> 


Scramjet lab 


Optics separate and 
Optics focus and collimate signal 
cross beam Flame/ 
Facility 


Nd:YAG laser (green) 
Broad-band dye laser (red) 
Narrow-band dye laser (ye!/ow) 


Beam Relay System 
Translates the measurement point 
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Setup at UVa: Optical System 2 au 


F CARS windows 
CARS | CARS windows 


[seratbeam | 


” oo 


From Laser Cart 


Dual-mode Scramjet 
Collection 
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@ Facility / Dual Mode Scramjet at UVa Sa 


Mach 2 nozzle Fuel injector CARS messpremen plane 
I | I I I I 
I y 


Nozz. & 


o 1 2 #3 4 


lsolator Combustor Constant Extender 
Facility z Area 


Air = 


Heater ft 
——se=_ nce gy) 


CARS measurement plane 


° Electrically heated clean air, nominal T,= 1200 K, MW = 2 nozzle 

¢ Long test times allow large data sets at steady inflow conditions 
¢ Highly turbulent flowfield 

¢ Good optical access 
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@) Nonequilibrium: CARS 


Mean 
Rotational 
Temperature 


=0.18 


Trotav: 


_s 


r 


rot 


and T,, 


600 800 1000 1200 1400 1600 1800 2000 


eS 


Mean 
Vibrational 
Temperature 


~=0.18 


GW 


UNIVERSITY 
WASHINGTON DC 


H, fuel, Config. C 


A. D. Cutler, G. 
Magnotti, L. Cantu, E. 
Gallo, P. M. Danehy, R. 
Rockwell, C. Goyne, 
and J. McDaniel, “Dual- 
Pump CARS 
Measurements in the 
University of Virginia's 
Dual-Mode Scranjet: 
Configuration C,” Paper 
AIAA-2013-0335, 2013. 


T 2vip OOServed to be 600 K higher than T,,, in facility freestream; not anticipated before experiment 


Vibrational Temperature 
CARS vibrational O, and Np 


— Vulcan, Tr — Vulcan, Tv 
-—#—CARS, Tr =<—CARS, TvN2 
—CARS, TvO2 -—®-CARS, Tv 


temperatures “lumped” 
together and used to 
validate the vibrational 
nonequilibrium model 
implemented in VULCAN 


0.0 
y/h 


-1.0 -0.5 0.5 


code by R. Baurle (NASA). 


A.D. Cutler, L. M. L. 
Cantu, E. C. A. Gallo, 
R. Baurle, P. M. 
Danehy, R. Rockwell, 
C. Goyne, and J. 
McDaniel, 
“Nonequilibrium 
Supersonic 
Freestream Studied 
Using Coherent Anti- 
Stokes Raman 
Spectroscopy,” AIAA 
Journal (in press 50 


@ Temperature Mean, Standard Deviation a4 


, Mean 
Rotational 
Temperature 

~=0.18 

~=0.49 

J . Trotav: 600 800 1000 1200 1400 1600 1800 2000 
T,or higher in free stream for p=0.49 due to ram mode 
Standard deviation 
~=0.18 
~=0.49 


El he 


% | ae sdTrot: 20 60 100 140 180 220 260 300 340 380 420 


Contours of Mole Fraction N, 


e Since inert shows overall mixing 


Mach 2 nozzle Fuel injector CARS i plane 
I | I I I I 
I I y 


Nozz. © 


0 1 2 3 4 


View from right 
looking 
upstream 


~=0.18 


=0.49 


N2av: 0.10.2 0.30.40.50.6 0.7 0.8 
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@) CARS compared with CFD (LES from NC State) GW 


Figure 9. Comparison between CARS and O, concentrations and RMS fluctuations for ‘0.17’ equivalence ratio. 


¢ CARS provides valuable quantitative comparison data 
— Fulton, Edwards, et al, (Configuration C) AIAA-2013-0117, 2013 


— The comparisons are not always this good! 
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6 CARS for nonthermal plasma AM 
¢ 100 Torr pulsed discharge at Ohio State 


v=0 Relative population 
| 0 
4 | 10 
@ @ @ i usec 
— 100 nsec |: = aa 
| 2 psec 
0.8} 1 usec Ls ss © © © 10ps60 
==""10 psec |: 10° — 
> 7 
0.6 sus udetsebeen' snvestseucesdncsesvecbectesdwsdesebcescesusttecceerecatesueodlevescea med ess 
= 
5.04 107 
” 
cre eg ee 
10° 
0 fevemefthanandicne ions See 
2100 2200 2300 0 2 4 6 8 10 


Vibrational quantum number 


(a) Raman Shift [cm] (b) 


¢ CARS can measure individual state populations for 


model validation. 
A. Montello, Z. Yin, D. Burnette, I.V. Adamovich, and W.R Lempert, “Picosecond 


CARS Measurements of Nitrogen Vibrational Loading and Rotational/Translational 
Temperature in Nonequilibrium Discharges,” J. Phys. D: Appl. Phys., Sept, 2013. 
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@ Additional Thoughts on CARS (a4 Gy 


CARS is a relatively accurate and precise high- 
temperature and species measurement technique in 
reacting and nonequilibrium flow 

— Temp. Accuracy 2-3%; Precision 2-3% (Magnotti, 2012) 
CARS can be performed through small windows or 
slots but need access on 2 sides. 

It is a relatively complicated technique, theoretically 
and experimentally, though seedless. 

— Takes months to analyze data carefully. 

ls being extended to much faster measurements ~10 
kHz with fs CARS 

— Very recently: 2 beam CARS, CARS imaging (Sandia) 

Has been performed using optical fibers to deliver the 
laser beams, both in ns and ps regimes ee 


| Conclusion S& 


¢ Variety of different molecular-based measurement 
technologies are available for measuring hypersonic 
nonequilibrium flows. 
— Can measure temperature, density, velocity, species 
concentrations, etc. 
— Which is best? Which technique to use? 
¢ Must start measurement campaign with an interview 
of “customer” to determine the requirements 
— How do different measurement technique capabilities 
satisfy the requirements? 
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o Conclusion S 


¢ Variety of different molecular-based measurement 
technologies are available for measuring hypersonic 
nonequilibrium flows. 
— Can measure temperature, density, velocity, species 
concentrations, etc. 
— Which is best? Which technique to use? 


e Must start measurement campaign with an interview 


of “customer” to determine the requirements 
— How do different measurement technique capabilities 
Satisfy the requirements? 
— Also consider practical issues: 
e Experience of the research team 
¢ Available equipment / budget 
¢ Schedule, difficulty of method, etc. 
— Rarely a perfect match, but a match can provide the data 
needed to advance or validate theory or computations 257 


| The interview: prior to getting started S 


¢ What parameter(s) need to be measured? 

¢ Must multiple parameters be obtained simultaneously to determine 
correlations? 

¢ What spatial resolution is required? 

¢ Is imaging required or are single-point or line measurements sufficient? 

¢ What temporal resolution is required? (e.g. time required for a single 
measurement) 

¢ Do measurements need to be time resolved? (e.g. a continuous sequence 
of data) 

e What accuracy is needed? 

e What precision is needed? 

¢ What quantity of data is required? 

¢ When is the data needed? Is instant (real-time) data required? 

¢ Where in the flow are measurements required? (inflow, exit, near walls, 
etc.) 

¢ What type of optical access is available? 

¢ Can (toxic) seed gases be introduced? Will they influence the properties 
being measured? 

¢ What is the ordered priority of the above requirements? 258 


